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ABSTRACT 

We present population synthesis calculations of the Tidal Downsizing (TD) hypothesis 
for planet formation. Our models address the following observations: (i) most abundant 
planets being Super Earths; (ii) cores more massive than ^ 5 — ISMq are enveloped 
by massive atmospheres; (iii) the frequency of occurrence of close-in gas giant planets 
correlates strongly with metallicity of the host star; (iv) no such correlation is found 
for sub-Neptune planets; (v) presence of massive cores in giant planets; (vi) gas giant 
planets are over-abundant in metals compared to their host stars; (vii) this over¬ 
abundance decreases with planet’s mass; (viii) a deep valley in the planet mass function 
between masses of ^ 10 — 2 OM 0 and ^ IOOM 0 . A number of observational predictions 
distinguish the model from Core Accretion: (a) composition of the massive cores is 
always dominated by rocks not ices; (b) the core mass function is smooth with no 
minimum at 3M0 and has no ice-dominated cores; (c) gas giants beyond 10 AU are 
insensitive to the host star metallicity; (d) objects more massive than ^ lOMjup do 
not correlate or even anti-correlate with metallicity. The latter prediction is consistent 
with observations of low mass stellar companions. TD can also explain formation of 
planets in close binary systems. TD model is a viable alternative to the Core Accretion 
scenario in explaining many features of the observed population of exoplanets. 


1 INTRODUCTION 


Core Accret i on m odel (CA; e.g., IPolIack et al.l 1 19961 : 
lAlibert et all l2005l ) stipulates that all planets grow from 
planetesimals - rocky or icy bodies ~ 1 km or more in size 
dSafronovI Il972l ). Planetesimals combine into bigger solid 
bodies by sticking collisions. More recent work suggests 
that planetesimals may have formed via streaming ins tabil- 
ities (lYoudin fc GoodmarJ l2005l : I Johansen et all l2007l ) and 
were born big, e . g., as large as 100 to 1000 km in size 
llMorbidelli et al] 20091 ) . In addition to this, pebbles, which 
are grains that have grown to the size of 1 mm to a few cm, 
are now s uspected to cont ri bute t o the growth of the cores 
strongly (lOrmel fc Klahi _ 2010|;_ Lambrechts fc JohansenI 


strongly llUrmel &: Klahil IzUlU: ILambrechts 
I2OI2I : [chamber J 20141 : Lambrechts et al. 20141 ). 


Whatever the growth mechanism of the cores, those 
that become massive attract gaseous atmospheres from pro¬ 
toplanetary discs. The atmosphere eventually becomes as 
massive as the core when the core mass exceeds a critical 
value. Merit ~ 10 M®, although the exact critical core mass 
is a function of dust o pacity, planet’s location and other im¬ 
portant physics (e.g.. fPerri fc Cameronlll974l : iMizun 3 ll98d: 
it evenso 3 119821 : lRafikovll2006l) . At this point a runaway ac¬ 
cretion of gas onto the core takes place, forming a gas giant 
planet (iPollack et al.l[l996l : iHubickvi et al.ll200a) . 

CA is the most widely accepte d theory of planet for¬ 
mation (e.g., see lHelled et al.ll2013l ). CA popularity is mo¬ 
tivated by its successes and, in no small measure, by the 


failures of th e alternative model. Gravitational disc Instabi l- 
ity (GI; e.g.. lKuipei[ll95ll : ICameron et al.iri982l : IBossII 19971 ). 
In particular, the classical version of GI cannot account for 
(i) the existence of terrestrial/rocky planets; (ii) any plan¬ 
ets within the inner ~ tens of AU of the host star; (iii) 
presence of massive cores inside of and metal overabun¬ 
dance of gas giant planets, and (iv) a positive giant planet 
frequency of occurrence - host star metallicity correlation 
dFischer fc Valenti l2005l ). CA, in contrast, has features (i- 
iii) built in by construction and predicts (iv) naturally as a 
result of producing more massive cores in high metallicity 
environments, s o that the runawa y gas accretion phase co m- 
mences earlier jlda fc Li3l2004c] f3: ]Mordasini et al.ll2014l ). 

However, a new planet formation framework called 
Tidal Downsizing (TD), in which GI is only the first step , 
has been suggest ed relatively recently llBolev et al.l I2OIOI : 
lNavakshirJl2010a ). The theory is an offspring of the Grav- 
itational disc Instab i lity model for planet formation (e.g., 
I Cameron et al. j[T13;[B OS 3 ll997l) but is far richer in terms 
of physics included in it. In a way, TD theory is GI the¬ 
ory modernised by physical processes many of which became 
standard features of GA at various times, but were somehow 
forgotten to be included in GI. These processes are planet 
migration, solids coagulating into massive bodies, pebble ac¬ 
cretion, and fragment disruption when compromised by too 
strong tidal forces from the host star. For a recent review 
of issues surrounding TD and GI, see iHelled et al.l (l2013l) . 
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and section 2 in iNavakshinI (l2015dl '). the latter specifically 
focused on TD. 

Several ways of addressing problems (i-iii) in the con¬ 
text of TD were qualitatively clear from its inception, but 
it is only very recently that a process accounting for (iv) 
was found. We have recently developed detailed population 
synthesis models of the new scenario enabling quantitative 
comparisons to observations. This paper, third in a series, 
presents a number of such comparisons and makes observa¬ 
tional predictions that may discriminate between CA and 
TD in the future. 

The origins of th is new theory lie in the apparently for¬ 
gotten suggestion of iKuin^ (Il95lli that GI may form not 
only the gas giant planets of the Solar System but also the 
rocky ones. He suggested that the inner Solar System plan¬ 
ets are made by destroying ~ Jupiter mass gas fragments 
within which du st sedimented down dMcCrea fc WilliamsI 
1 19651: Ib osslll998h to form massive and dense cores composed 
of heavy elements. Hydrogen/helium and other volatile com¬ 
ponents of the fragments are disrupted by the Solar tides and 
eventually consumed by the Sun, whereas the much denser 
cores survi ve to becom e the p resent day planets. 

Until iBolev et al.l ll20ld l. no physical way of actu¬ 
ally placing the massive gas fragments at ~ a few AU 
dista n ces from the Sun seemed to exist (e.g., I Rice et ^ 
I 2 OO 5 I: iRafikovI [ioosib This appears to be the main rea¬ 
son why this avenue of planet forma tion was discounted 
early on (iDonnison fc WilliamsI Il975l l. However, we now 
know that the fragments do not have to be born 
wher e the planets are now because of planet migratio n 
fe.g..lLin fc PapaloizorJfl979l: iGoldreich fc Tremainef 


iBolev et al.l 1 2010h found that ~ Jupiter mass gas fragments 
born by gravitational instability in the outer {R ~ 100 AU) 
cold protoplanetary disc do not stay there, as usually as¬ 
sumed, but migrate inward rapid l y (in as little as ^ 10'* 
yrs; see al so Baruteau et al.l 1201 ll : Icha fc NavakshinI I 2 OIII : 
IZhu et ^ 12012 ! '). The fragments are initially very fluffy , 
and it takes up to a few Myrs (e.g., BodenheimeJ Il974 
iBodenheimer et al.lll980l : IVazan fc Helledll2012l ~) for them to 
contract and collapse to what is usually taken as the time 
t = 0 configuration of G I planets (the so ca lled ’’hot start” 
of gas giants, see, e.g., iMarlev et al.l uOOn b Since migra¬ 
tion process is much quicker, the pre-collapse gas fragments 
are usually disrupted by tides from their host stars when 
the fra gment separation shrinks to a few AU. INavakshinI 
ll2010al ~l arrived at similar ideas via analytical estimates and 
simulations of core form ation within isolated gas fragments 
llNavakshirJl2010bl . l20lih . 

Since then, several dozen investigations into the physics 
of TD or into processes important to TD (cf. f urther refer¬ 
ences in iHelled et al.ll2013 : INavakshinI l2015dh sprung up. 
There was however much less work on trying to relate the 
TD hypothesis to the observations of exoplanets in a quanti¬ 
tative way, which is not particularly s urprising given the the - 
ory is still in its infancy . Neve rtheless, iForgan fc Ricd (l2013l l 
and iGalvagni fc Maverl (l2014l l built the first population syn¬ 
thesis models for TD in order to compare the predictions of 
the models with the statistics of observed exoplanets (fo r 
the latter, see a recent review by I Winn fc Fab rvckvl 120141 . 
The results of these two studies are somewhat divergent. 
iForgan fc Ricd (l2013f) find that TD is incapable of producing 
any planets at small separations (R iS 10 AU) since most of 


their gas fragments are tidally disrupted before massive solid 
cores could form within them. iGalvagni fc Ma\^ (l2014l ~) did 
not include grain sedimentation in their models, hence could 
not say anything about the post-disruption core-dominated 
planets. However, their study found a wealth of Jovian mass 
planets that successfully migrated into the inner disc, avoid¬ 
ing tidal disruptions, and hence they suggested that TD may 
well be effective in producing a hot-Jupiter like population 
of gas giants. 

Another challenge to TD is the expectation that higher 
metallicity fragments are more likel y to be tidally disrupted 
becau se of slower radiative cooling (iHelled fc Bodenheimerl 
I 2 OIIII . This would imply that low metallicity environments 
must be more hospitable to gas giant formation via TD chan¬ 
nel. Since a strong positive giant planet frequency - metallic¬ 
ity co rrelation is observed (lGonzalejll999l : [Fischer fc Valentil 
l2005h . the general feeling is that TD is strongly disfavoured 
by the data. 

A solution to these and other sh ortcomings of TD 
hypothesis may be ’’pebble accretion” llNavakshinIl2015al l . 
the process in which large 1 mm to a few cm) 

grains from the disc separate out of the gas-dust flow 
past an embedded massive bod y (large planetesimal or a 
planet) and accrete onto it (e.g., Johansen_j£Lacerdal l201(]l : 
lOrmel fc KlahrlI 2 OIOI : iLambrechts fc Johanserj 2012ll . Peb¬ 
ble accretion is mainly called upon in the Gore Accretion 
model to ac celerate the growth of solid cores at large dis¬ 
tances (e.g.. IHe lled fc Bodenheimerl 1201^ : IChambera l2014l : 
ILambrechts et ah 2014h : no application of pebble accretion 
to TD was however done until very recently. 

A ID radiative hydrodynamics study with grains 
treated as a second fluid showed that external pebble depo¬ 
sition in pre-collapse gas fragm ents significantly a ccelerates 
their contraction and collapse (INavakshinI l2015'r] L Using a 
ID viscous disc evolution code to treat the disc-planet in¬ 
teractions (angular momentum and heat exchange, pebble 
accretion) and co-evolution, we confirmed that pebble accre¬ 
tion allows more fragments to survive tidal disruptions, and 
also does so prefer entially in hi gh metallicity env ironments 
llNavakshinll2015'9l . In paper I (INavakshinI 1 201 5^ . this nu¬ 
merical scheme was extended to include grain growth, sedi- 
menta tion and core form ation within the fragments. In pa¬ 
per H llNavakshinll20153) . a large set of population-synthesis 
like experiments in TD setting was performed. Since peb¬ 
ble accretion strongly increases grain and metal abundance 
within the gas fragments, massive (Mcore ~ 10 M 0 ) cores 
composed of heavy elements form much quicker than they 
do inside “pristine” fragments made of the dust/gas compo¬ 
sition equal to that of the host star. For this reason, an abun¬ 
dant massive core presence in the inner dis c was found in 
paper II, in contrast to the earlier results bv IForgan fc Ric'3 

Interestingly, while survival of gas giant planets in 
these numerical experiments was strongly enhanced at high 
metallicities, the same was not true for core dominated 
planets with mass iS 10 M 0 . These planets were most 
abundant around ~ Solar metallicity hosts rather than at 
higher metallicities. This result is in line with both ra- 
dial velocity and transit observations of exoplanets (e.g. , 
ISousa et al.ll2008l : iMavor et al.ll20jTl : iBuchhave et al.ll201^ 
which showed that planets smaller than Neptune in mass 
or radius do not correlate with metallicity of the host. 
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Strictly speaking^ this contradicts expectations from CA 
dlda fc LirJ l2004a i that massive cores are more abundant 
at high metallicities. This is the feature of the model called 
upon to explain the gas giant planet positive metallicity cor¬ 
relation with metal abundance of th e host stars. While more 
recent work l|Mordasini et al.ll2012l 'l shows that CA can ex¬ 
plain the absence of metallicity correlation for Neptune mass 
or smaller planets, this result depends sensitively on uncer¬ 
tain details of type I migration prescription (see section 7.2 
in the paper quoted just above). 


Here we improve on methods of paper I and II in a 
number of ways, inyestigate the results in greater detail, 
and present additional model predictions to distinguish TD 
from CA observationally. First, we study a much broader 
range of the initial fragment masses, e.g., from the minimum 


of MminO = 1/3 M j to the maximum of MmaxO = 16 Mj 
l|Navakshinl l2015d . studied fragments in a much narrower 
mass range, 0.5 < Mq/Mj < 2). This is important as the 
initial fragment mass cannot be accurately constrained at 
this stage: different authors arrive at d ifferent answers (e.g., 
Bolev et al.l 20ld : For gan fc Ricellioill '). Secondly, following 


Navakshin et al.r 1 2013 1. we calculate the structure of the 


dense gas layers adjacent to the core growing inside the gas 
fragment, hence resolving the gravitational influence of the 
core on the very centre of the gas fragment. This step is 
key to studying Neptune to ~ sub-Saturn mass planets. 
In the context of TD, such planets form when the core is 
surrounded by massive gas atmosphere bound to the core. 
This is physically similar to how massive ga s atmospheres 
form around cores in the cont ext of CA (e.s.. iMizun 3 I1980I: 
IStevens^ll982l : lRafikovl2006l l. but all the action takes place 
in the centre of the massive gas fragment (formed by GI and 
not due to the presence of the core) rather than the pro¬ 
toplanetary disc. When the gas fragment is disrupted, the 
bound atmosphere survives together with the core since it 
is far denser than the rest of the fragment ll Navakshin et all 
|2014 1. Finally, we now use a more realistic - better observa¬ 
tionally constrained - approach to the protoplanetary disc 
dispersal to end our runs. 


We find that our population synthesis model reproduces 
a number of observational facts/correlations for planetary 
properties that were previously suggested to vindicate CA 
as the only correct theory of planet formation, such as the 
metallicity correlations, the presence of cores inside of and 
the general overabundance of metals in gas giant planets 
compared to their host stars, the existence of a critical core 
mass above which the core must possess a massive envelope 
composed of volatiles and hydrogen/helium. Fortunately, 
there turn out to be a number of observational predictions 
distinguishing our theory from CA model. 

Section 2 of this paper presents a summary of the main 
assumptions and numerical methods used here. Readers not 
interested in computational detail of the model may skip 
section 2 and proceed to §3, where a broad brush overview 
of the results is presented. Planet mass function (PMF), 
metallicity correlations, planet and core compositions are 
described in §4-9. The radial distribution of simulated plan¬ 
ets is compared to observational constrains in ESI A broad 
discussion of the implications of the results is given in m 


2 COMPUTATIONAL APPROACH 
2.1 Main assumptions 

Our main assumptions are very similar to those in paper I 
and II. We recount them briefly below, with more detail on 
the calculations to follow in 112.21 below. 

(i) We study the migration of a gas fragment in an ini¬ 
tially massive self-gravitating disc around a star with mass 
M* = 1 Mq in the last phase of the disc existence. The 
fragments are born at separation a ~ 100 AU where this 
process is allowed by physical conditions. During our sim¬ 
ulations, the mass supply from the parent molecular cloud 
has ceased, and the disc mass monotonically decreases due 
to accretion on the star and the disc photo-evaporation. This 
approach neglects the earlier population of fragments born 
when the star was less massive for the sake of simplicity at 
this exploratory stage. 

(ii) Only one fragment per disc is simulated. This is a se¬ 
rious shortcoming of the calculations as simulations of self- 
gravitating discs show formation of multip le fragments form¬ 
ing and interacting with each other (e.g., lBolev et al.ll20ig : 
ICha fc NavakshirjT201ll : iMerul [201^ . We attempt to miti¬ 
gate for this by variations in the initial parameters of the 
disc and in the fragment’s migration speed (see 12.21) . 

(iii) Planetesimals are neglected in this paper. These mi¬ 
nor bodies are a side show in the context of TD, something 
that also happens during planet formation but is not a pre¬ 
requisite to formation of any of the planets. Moreover, in TD 
planetesimals are formed by the proto-planets (that is, mas¬ 
sive g as fragment s ) and not the other way around. Specif¬ 
ically, iNavakshinl (1201 il l argued (in §7 of that paper) that 
environments inside massive gas fragments are both bet¬ 
ter and safer for solids to grow than the main body of the 
protoplanetary discs. Higher gas densities and self-gravity 
of the gas fragments provide high grain growth rates and 
protection against fragmenting high speed collisions, turbu¬ 
lence and the rapid migration of the fragments into the sta r 
(the famous ”1 metre b arrier”, see IWeidenschillind Il980l l . 
iNavakshin fc Ch^ (l2012l l showed that within the fragment, 
solid bodies with small angular momentum accrete onto the 
central massive core, whereas ^ 1 km or larger bodies with a 
larger angular momentum orbit the core instead. When the 
gas fragment is disrupted, bodies closest to the core remain 
bound to it, possibly contributing to formation of satellites 
of future super Earth or more massive planets. Bodies far¬ 
ther away from the core become gravitationally unbound 
from it when the fragment is disrupted. They however re¬ 
main bound to the parent star and form debris rings with or- 
bital properties not unlik e the Kuiper and the asteroid belts 
jNavakshin fc ChcJl2012l ~). Most importantly, there is likely 
much less mass in these fragment-produced ” planetesimals” 
than in the Core Accretion scenarios of planet formation. For 
this reason, although planetesimals formed in earlier gener¬ 
ations of disrupted fragments may accr ete onto and influ¬ 
ence the gas fragm ent evolution (e.g., see lHelled et al.ll2008l : 
iBolev et al.ll201ll i. this process is neglected in the present 
paper. 

(iv) Once formed by gravitational ins tability, fragments 
do NO T accrete more gas from the disc. INavakshin fc Cha 
1I2OI3II performed 3D SPH simulations of gas fragments em¬ 
bedded in massive marginally stable discs that included 
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a physically motivated analytical prescription for radiative 
heating of the gas around the planet’s location. It was found 
that planets less massive than ~ 5 — 10 Jupiter masses cre¬ 
ate a hot atmosphere around themselves that stifles accre¬ 
tion of more gas onto them. More massive planets accrete 
gas rapidly and run away well into the BD regime at which 
point migration of the fragment is stalled. More 3D simula¬ 
tions are needed to ascertain (in)dependence of this result on 
dust opacity of the disc and other parameters of the model. 

(v) Pebbles are accreting onto the fragments at the rate 
(equ ation 5 in paper II) calculated by extrapolating results 
from iLambrechts fc JohansenI (l2012l l. The size of the pebbles 
is set to 0.3 mm for simplicity, and the incoming pebbles are 
deposited into the outermost layers of the planet. 

Readers not interested in technical details of the calcu¬ 
lations may skip section [2721 and continue directly to section 
[3 where results of this paper are beginning to be presented. 


2.2 Summary of numerical methods 

We follow numerical methods explained at length in paper I 
and II. These are summarised below. We also explain several 
important modifications that improve our methods. 


2.2.1 Disc evolution 


We set the inner boundary of the disc to Rm = 0.08 AU, 
whereas the outer boundary is set to Rout = 400 AU for 
this paper. We use 300 radial bins, logarithmically spaced 
in , which gives us the highest numerical resolution near 
the inner boundary. We found that increasing the number 
of radial bins to over 1000 only led to changes in the planet 
migration times at the level of 10%. We hence chose 300 
bins as acceptably accurate for a population synthesis study: 
many parameters of the model (such as viscosity parameter 
a) are known with a far worse precision anyway. 

The disc is init ialised with a surface density profile 
E(R) oc (1/R)(1 - yjRin/R) exp(-R/Rdisc), where Raise = 
100 AU is the disc radial length scale, with the normalisa¬ 
tion set by the total disc mass, Ma, integrated from Rin to 
Rout ■ The disc surface density is evolved by solving the time 
dependent viscous disc equation (cf. equation 1 in paper I), 
but now also including the disc photo-evaporation term to 
be explained below. The disc exchanges angular momentum 
with the planet via a prescription switching smoothly from 
type I into type II regime. The latter regime, when a gap in 
the disc is opened by the planet, is modelled by the widely 
used ’’impulse approximat ion” (e.g.. lLin fc PapaloizoulflOSB : 
lArmitage fc Bonnellll200^ i . The type I regime (no gap in th e 
disc) is motivated by the results of iBaruteau et al] ll201lli . 
parameterised via the planet migration timescale, tj, given 
by 


h = frr 


Mi 


( 1 ) 


where a is the current position of the planet, Md = vrEa^ 
is approximately the local disc mass, Qa is the local Kep- 
lerian angular velocity, and /migr ~ 1 ( see below) is a free 
parameter. The criterion introduced bv ICrida et al.l (l2006l i 
is used to switch the planet migration between type I and 
type II. We emphasise that our approach conserves angular 


momentum explicitly: the torque from the disc on the planet 
is equal to that of the planet on the disc with the minus sign, 

of course. _ 

_ The disc viscosity is modelled as a lShakura fc SunvaevI 

lll973tl prescription, with the a viscosity parameter given by 
the sum of a constant term, ao = 0.005 and a term that 
depends on the local self-gravity of the disc. 


Qsg 


= 0.2 


Qo 

Ql + Q^ ' 


( 2 ) 


where Q is the local Toomre’s parameter. The parame¬ 
ter Qo, setting the tran sition from a viscosi ty driven by 
gravito-turbulence (e.g., iLodato fc Ricd l2005l i to the non 
self-gravitating one, is set to Qo = 5 here. We found that 
lower values of Qo (Qo = 2 was used in paper II) result in 
the stars with discs fraction that is too high at early times, 
t ^2 Myrs. To fit observational constraints, one would ei¬ 
ther have to choose a much larger value of ao or require 
implausibly high photo-evaporation rates. 

In paper II, disc photo-evaporation due to the host 
star’s UV and X-ray radiation was modelled as the sum 
of th e fits to the respective photo- evap oration rates Sin(R ) 
from [Alexander fc Armitagef (l2007f) and I Owen et al ]iS. 
The UV ionising photon luminosity of the star was set to 
"hion = photons s“^, while the X-ray flux of the star was 
Lx = 2 X 10^° erg s“^. In this paper we employ the same ap¬ 
proach, but also include the external photo-evaporation rate, 
Eext, by UV flux from nearby mass ive stars at 0.1 times the 
rates parametrised bv lClarl3 (l2007l i. This procedure in gen¬ 
eral is not guaranteed to satisfy ob servational const r aints o n 
disc dispersal, however. Following iMordasini et al] (l2009al l. 
we introduce a free parameter, 0 < Cev, so that the photo¬ 
evaporation rate is given by 


Eev — Cev ^Ein -f Eext^ 


(3) 


With this photo-evaporation prescription, the total photo¬ 
evaporation rate fr om the disc is equal to ARv = 3.1 ^ev x 
10 “® Mq yr“^. As IMordasini et ^ ll2009ah . we generate a 
uniform random distribution of log ^ev bounded from below 
by Cmin and ((max, and evolve a sample of ~ 1000 discs with¬ 
out planets. We then plot the fraction of stars with discs 
as a function of time, and compare it with the observation- 
ally e stablished depende nce, ~ exp{—t/tob), where tob = 3 
Myr (iHaisch et al.ll200ll '). By experimenting with values of 
Cmin and Cmax, it was found that a reasonably good agree¬ 
ment with the observations is obtained for ((min = 0.02 and 
Cmax = 10 (cf. hg. [T]). There is a certain degeneracy in this 
procedure. A higher value of disc viscosity (ao) could be off¬ 
set by reducing the photo-evaporation rate, but our main 
results do not depend on this level of detail. 


2.2.2 Planet evolution module 

The planet evolution module used in this paper is nearly 
identical to that described in section 4 of paper I. In brief, 
we model the pre-collapse gas fragment contraction using 
the “follow the adiabats” approach, which assumes the frag¬ 
ment to be strongly convective except for the boundary, 
and includes the radiative cooling of the fragment mod¬ 
ified by the irradiation it receives from its surroundings 
(which can be very important, even unbinding the planet. 

















































Tidal Downsizing III 5 


Disc Evolution, a = 0.005 



Time [Myr] 

Figure 1. Fraction of stars with discs as a function of time for 
several values of the minimum disc evaporation factor fmim as 
labelled in the legend. The solid line gives the curve exp(—t/iot,), 
where tab = 3 Myr. Cmin = 0.02 provides an acceptable fit to the 
solid curve. 


see ICameron et al.l Il982l : IVazan fc HelledI l2012l ^ . The frag¬ 
ment is embedded in the disc if it migrates in the type I 
regime, and is hence bathed in the disc radiation with effec¬ 
tive temperature of the disc mid plane. When the fragment 
migrates in type II regime, the irradiation from the central 
star and the edges of the gap carved by the fragment may 
be important. One difference from paper I is that we modify 
the intrinsic luminosity of the planet, Trad, by introducing a 
“green house” effect for the planet immersed into the disc: 

L,ad = ^ , (4) 

where Liso is the luminosity of the fragment in isolation, 
and Td = kEci/ 2 is the disc optical depth at the planet’s lo¬ 
cation. This correction is necessary because the irradiation 
luminosity incident on the planet from the disc was calcu¬ 
lated in paper I and II neglecting the planet’s energy input 
to the disc around the planet’s location. Furthermore, the 
disc thermodynamic properties are azimuthally averaged in 
our ID approach, which additionally washes out the planet’s 
influence on the thermal balance of the disc in the planet’s 
vicinity (~ Rh or H, whichever is larger). The planet’s in¬ 
put into the disc temperature around its location must how¬ 
ever be taken into account since massive planets (several 
Jupiter masses or more) can easily dom inate the disc en¬ 
ergy production around the planet le.g.. iNavakshin &: Chal 
l2013l f . The planet’s luminosity calculated via equation |4] is 
therefore used in place of Liso in equation (9) of paper I to 
determine the rate of planet’s radiative energy loss. 

The dust physics used in this paper is identical to that 
described in sections 4.2, 4.3, and 4.6 of paper I. We use 


three grain species - water, CHOnQ and rocks in relative 
abundances of 0.5, 0.25 and 0.25. 

The total grain mass abundance is initially equal to 
the metallicity of the system, Z. The grains grow by stick¬ 
ing collisions, can sediment down modulo turbulence and 
convective grain mixing which oppose grain sedimentation. 
Grains can also be fragmented if they sediment at too high 
velocity, which in practice limits their size to a few cm at 
most. Grains vaporise once the given species vaporisation 
temperature is exceeded. Gore assembly by grain sedimen¬ 
tation therefore depends sensitively on the temperature and 
convection of the inner regions of the gas fragment. The 
core growth prescriptions are identical to those presented in 
section 4.4 of paper 1. 


2.2.3 Atmospheres bound to solid cores 

One significant modification to our treatment of the planet’s 
internal evolution is addition of an “atm osphere” around 
the gr owing massive cores. As found by INavakshin et al.1 
(I2OI4II . when the core reaches the mass of a few to a few 
tens of M0, the weight of the gas surrounding the core 
becomes comparable t o the weight of the core itself. In 
INavakshin et akl (l2014ll . the structure of the core’s atmo¬ 
sphere was calculated assuming hydrostatic and thermal bal¬ 
ance. Additionally, the outer boundary conditions for the 
atmospheric structure calculations were fixed by joining the 
atmosphere to a non-perturbe d approximate ana lytical so¬ 
lution for a gas fragment from iNavakshinI (l2010al l. Unpub¬ 
lished experiments done with a ID radiation hydrodynam - 
ics code (described most recently in lNavakshinll20i3 . l2015^ 
showed that collapse of the atmosphere around the core does 
not actuall y lead to the collapse o f the whole fragment, as 
assumed in INavakshin et ^ (l2014f l. The physical reason for 
this is that during such a collapse contraction of the atmo¬ 
sphere generates a significant compressional energy, which is 
then communicated to the rest of the fragment by radiation 
and convection. The fragment then expands, which lowers 
the density at the fragment-atmosphere boundary. As the 
result, the core’s atmosphere’s outer layers also expand out¬ 
wards, preventing further collapse of the gas layers near the 
core. This “feedback loop” prevents a hydrodynamical col- 
lapse of the who l e frag ment that was suggested to occur in 
INavakshin et al.1 (12014 1. Physically, we expect a slow con¬ 
traction of the atmosphere instead at the rate regulated by 
the cooling of the whole fragment. 

To take this physics into account, we calculate the mass 
(and structure) of the atmosphere near the core at regular 
(short) time intervals. The atmosphere region is defined as 
that around the core where the total specific gas energy, 
u — GM{r)/r < 0, where u is the specific internal gas en¬ 
ergy, r is the radius from the core, and M(r) is the enclosed 
mass within that radius, including both gas and the core. If 


^ CHON is a mnemonic acronym for a material dominated by car¬ 
bon, hydrogen, oxygen, and nitrogen, excluding water. Its com¬ 
position, material properties and vapor pressure are here taken 
to b e similar to t hat of the grains in the coma of Comet Hal¬ 
ley ‘ job erd ll2004r) . CH ON is a widely used component in pla net 
formation theory fe.g.. IPollack et ^ll996ldHelled et al.ll2008l) . 
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the fragment collapses due to H 2 dissociation then no mod¬ 
ifications are done, since the atmosphere just becomes the 
part of the gaseous envelope of the planet. However, if the 
fragment is tidally disrupted, then the mass of the disruption 
remnant should be more massive than just the core (which 
was assumed in papers I and II). Specifically, if at disruption 
the mass of the atmosphere, Matm, is lower that the critical 
mass for atmosphere collapse, Mac = (l/3)Mcore, then the 
mass of the survived planet is given by Mp = Mcore -b Matm- 
However, if Matm > Mac, then we expect that a fraction of 
the total mass of the gas envelope would have accreted onto 
the core. Since our planet formation module is not yet so¬ 
phisticated enough to calculate this fraction from first prin¬ 
ciples, we simply endow the core with gas mass randomly 
distributed in log between Mac and the total mass of the 
fragment minus that of the core. This procedure is equivalent 
to a partial tidal disruption of a fragment, which typically 
results in planets with masses between a few tens of Mq to 
~ 1 Jupiter mass. We find that partial tidal disruptions of 
fragments more massive than ~ 2 Mj are rare. 


2 . 2.4 Monte Carlo initial conditions 

Table 1 summarises the randomly drawn variables and their 
ranges used in this paper. The approach is similar to that 
from paper H: for each random variable, the distribution is 
uniform in the log of the parameter, and is distributed be¬ 
tween the minimum and the maximum values given in Table 
1. For example, the initial mass of the gas fragments is Mq in 
Jupiter masses, with the log Mq distributed randomly from 
log(l/3) to log 16. Note that the logarithmic distribution for 
Mq translates into the dN/dMo oc I/Mq dependence for the 
mass function of the fragments, which corresponds w ell to 
the one observed for gas giants (|lJdrv fc Santosll2007l . finds 
dN/dMo oc M|j"^ °®). This range is far broader than the one 
used in paper H, where we only sampled the planets with 
initial masses in the range 0.5 Mj < Mq < 2 Mj. We hope to 
make tentative connections of TD hypothesis to the brown 
dwarf regime. 

The disc mass is sampled between 0.075 Mq and 
0.15 Mq of that. As explained in paper II, the mass of the 
disc at which the disc is gravitationally unstable (Toomre’s 
parameter Q = 1.5) is about 0.15 Mq (cf. fig. 1 in paper 
H) for the chosen initial surface density profile. Lower disc 
masses are considered here because 3D simulations show it 
is possible for fragments to be born in less massive non self- 
gravitating discs due to fragmentation induced by pertur¬ 
bations fro m fragments born at previous more massive disc 
epoch (e.g.. lMerull2013l L The fragments initial locations are 
sampled between 70 AU and 1.5 times that, which is reason¬ 
able since the minimum in the Toomre parameter is reached 
at a ~ 90 for our disc. 

The opacity reduction factor is se t to 1 in this p aper 
(that is, the interstellar grain opacity of IZhu et al.ll200^ . not 
modified by grain growth, is used). The planet migration 
factor, /migr, is sampled between 0.5 and 2, which is much 
smaller a range than the 1 to 10 range used in paper H. 
As in previous papers, we do not actually model the inter¬ 
nal structure of the core, assuming it has a fixed density of 
5 g cm“®, and we hence cannot calculate the core’s lumi¬ 
nosity self-consistently. To model the core’s luminosity, the 
Kelvin-Helmholtz contraction time of the solid core, tkh, is 


defined isee lNavakshin et al.ll20l3 l. and the core’s luminos¬ 
ity is prescribed as in section 4.4 of paper I. The parameter 
tkh is varied between 10® and 10^ years for the simulations 
in this paper. 

The only exception from the uniform random param¬ 
eter distribution procedure is the host star/disc metallic- 
ity distribution, which can be constrained observationally. 
It appears to be close to a Gaussian dist ribution in the 
CORALIE plan e t samp le dUdrv et al.l 12003 1 . As shown by 
iMordasini et al.l (|2009al L a good fit to the observations is 
provided by the Gaussian distribution 


dp{ZL 

dZh 


r(27r)i/' 


■ exp 


2cr2 


(5) 


where Z^ =[Fe/H] = [M/H], the usual logarithmically de¬ 
fined metallicity, a = 0.22, and dp/dZr is the probability 
density. The subscript “L” on Zl stands for “logarithmic”, 
since Zl = login(-^/^Q), where Zq = 0.015 is the Solar 
metallicity llLoddersll2003ll . e.g., the fraction of mass in as- 
trophysical metals compared to the total mass. 

The metallicity distribution of the host stars used here 
is hence observationally motivated and differs from the hxed 
metallicity bins used in our previous papers. The metallicity 
of the disc and the initial metallicity of the planet, Zq are 
equal to that of the host star. 


Zq — Zq 10^^ . 


(6) 


3 RESULTS OVERVIEW 

We performed 20,000 fragment-disc evolution experiments 
for this paper. In this section we provide an overview of the 
results, with further detailed analysis to follow. 


3.1 The mass-separation plane 

Figure [2] presents the population of planets produced by 
the end of our runs in the planet mass - separation plane. 
The mass of the planets is shown in the units of Jupiter 
masses on the left vertical axis, and in the Earth mass units 
on the right. To improve clarity of the figure, we reduce 
the number of planets (individual symbols) shown in the 
hgure by randomly selecting a sub-sample of the planets. 
The regions to the left of the vertical dashed line (set at 
0.1 AU) and below the horizontal dotted line (set at 0.05 
Mj) show only l/20th of the sample since these regions are 
very strongly crowded with planets. The remaining region - 
above the horizontal line and to the right of the vertical line 
- shows 1/2 of the sample. 

The planets to the left of the vertical line are those 
that migrated through our inner boundary condition (oper¬ 
ationally, they reached very close to it, a = 0.09 AU). Our 
choice of the inner disc radius, Rin = 0.08 AU, is somewhat 
arbitrary. It is likely that most of these planets will be driven 
so close to the star that they are either consumed by the star 
or are unbound by over-heating due to the intense radiation 
held of the star. However, we expect that magnetospheric 
interactions will cutoff the disc at different distances from 
the star in different systems. This cutoff distance will also 
vary with time in a given system as the accretion rate and 
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Table 1. The range of the Monte Carlo parameters of the population synthesis calculation. The first row gives parameter names, the 
next two their minimum and maximum values. The columns are: Planets initial mass, Mq, in Jupiter masses; (^ev, the evaporation rate 
factor (cf. eq.[3l); /p, the pebble mass fraction determining the fraction of the disc grain mass in the pebbles; ao [AU], the initial position 
of the fragment; the initial mass of the disc, in units of Mq] /migr^ fhe type I planet migration factor; Million years, determines 
the luminosity of the core; ad, turbulence parameter within the fragment; the grain breaking velocity, in m/s 

Parameter Mq Cev fp ao Qd 

Min 1/3 0.02 0.04 70 0.075 0.5 10® lO”"* 5 

Max 16 10.0 0.08 105 0.15 2 lO'^ 0.01 15 


the magnetic field intensity vary. Due to these effects, we ex¬ 
pect some of the planets that fell trough our inner boundary 
to stop somewhere in the region between the star and Rin- 
For this reason we shall include a fraction of these planets in 
our further analysis and hence we show them in fig. [2] For 
these planets, the separation is a uniform random variable 
distributed between 0.03 AU to 0.09 AU. 

The colours of the symbols are used to indicate the 
metallicity of the parent star as noted in the legend. For 
example, green symbols show the most metal poor stars, 
[M/H] < —0.25, whereas red show the most metal rich ones, 
[M/H] > 0.25. One can notice the prevalence of red over 
green for gas giant planets in the inner few AU , signifying 
the po sitive metallicity correlation reported in iNavakshirJ 
(l2015lJ) and paper II. The outer tens of AU gas giants 
and the sub-Neptune planets below the horizontal line show 
much less, if any, metallicity preference, hinting on the ab¬ 
sence of metallicity correlation there. These results were also 
found in paper II. 

3.2 Statistics of fragment destruction and survival 

We shall now look at the results in a different way. All 20,000 
of our fragments are initially born at 70 < a < 105 AU, yet 
there are just some hundreds of planets left in the region 
a > 10 AU by the end of the simulation. It is therefore clear 
from fig. [2] that most of the fragments migrate very close to 
the star, and are either disrupted and end up as low mass 
planets below the horizontal line, or migrate through the 
inner boundary of the disc and appear to the left of the 
vertical line. 

For a more quantitative analysis of this. Figure [3] sum¬ 
marises the fate of our gas fragment migration experiments 
in two different ways. Fig. [3^ shows the initial fragment 
and the final planet mass spectrum. The former is shown 
with the black histogram, and is uniform in mass from 
Mmino = 1/3 Mj to Mmaxo = 16 Mj. The quantity shown 
in the histograms is dNp/dMpAMp, the number of planets 
found in the mass bin between Mp and Mp -|- AMp. Some 
deviations from a fiat line in the figure are due to random 
Monte Carlo fluctuations in the fragment mass initialisation 
procedure, and also MminO and MmaxO falling inside the his¬ 
togram bins rather than precisely on the bin edges. 

The coloured histograms show the final masses of the 
planets after TD ’’processing” of the fragments in the disc. 
These are shown for all planets independently of where 
they end up in terms of their final planet-host separation. 
The populations are further broken on the metal rich bin, 
[M/H] > 0 (red colour), the metal poor bin, [M/H] < 0 
(green) and the total (blue). 

We see that the distribution of final planet masses 


broadly divides into two groups: the disrupted one (frag¬ 
ments less massive than 1/3 Mj or ~ 100 Mq) and those 
that avoided the disruption. The former population is more 
populous than the latter, as we shall see. The most populous 
group of disrupted planets are those that have masses from 
sub-Earth mass to ~ 10 M®. We notice that sub-Earth plan¬ 
ets are more abundant in metal poor discs whereas super- 
Earths are more abundant at higher metallicities. This dis¬ 
tribution however contains planets at all separations from 
the star and needs to be cut into more specific radial bins 
to learn more, which will be done below. 

To digest the overall results of our simulations further, 
we find it useful to group the end outcome of the runs into: 

(i) Disrupted planets - those that were tidally disrupted 
anywhere in the disc between Rin and their birth position. 

(ii) “Assimilated” planets: those that migrated all the 
way to the inner boundary of our disc, i?in. Our calcula¬ 
tions are oversimplified at the inner boundary, as pointed 
out above. In all likelihood, most of these fragments would 
be pushed closer to and possibly even into the star and thus 
be accreted (assimilated) by the star, although a minority 
may survive when the disc dissipates away. 

(iii) “Cold giants” - those fragments that survived tidal 
disruption and matured into a giant planet at large (a > 
10 AU) separation from the star. 

(iv) “Warm giants” - those fragments that survived tidal 
disruption, and, despite the rapid inward migration, were 
able to stop outside the inner disc boundary but within a < 
10 AU. 

The frequency of these four outcomes are displayed in 
fig. [Sjo separately for the four metallicity bins. Note that 
tidal destruction of the fragments is the most common out¬ 
come 80% of the time) for Solar metallicity or less 
metal rich stars, but becomes progressively less common 
for the two metal rich bins, dropping to just over 20% for 
[M/H] > 0.25. This result is related to the posi tive metal¬ 
licity c orrelation for gas giant planets found by iNavakshinI 
ll2015bl ) and in paper H: the higher the metallicity, the fewer 
giants are disrupted, and hence more giants survive per star. 

Most of the fragments in the most metal-rich discs (red 
curve in fig. avoid tidal disruption. However, most of 
them end up migrating through the whole of the disc to i?in. 
As already suggested, a fraction of them is likely to survive 
inside Rin as ’’hot Jupiter”. Focusing now on the population 
of gas fragments that were not tidally disrupted and not mi¬ 
grated through all of the disc (the last two entries on the 
horizontal axis in the figure), we observe that most of these 
reside at large separations from the star, a > 10 AU. The 
gas giant planets in the inner few AU is hence the rarest gas 
giant population in our synthetic models. Finally, it is pos- 
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Figure 2. Simulated planets in the planet mass versus planet-host separation plane. The symbols and colours reflect the metallicity of 
the host stars, as detailed in the legend. For the sake of clarity, only l/20th fraction of planets is shown to the left of the vertical line 
and also below the horizontal line. 1/2 of all planets is shown in the right top quartile of the plane. 


(a) Planet disruption outcomes 



(b) Gas Fragment Fate 



Figure 3. (a) Comparison of the initial mass spectrum of fragments (black diagram) to the final distribution of planetary masses. The 
colours group the results by metallicity, with low metallicities defined as [M/H] < 0 and high defined as [M/H] > 0. (b) The outcome of 
gas fragment migration experiments in terms of the fraction of the initial fragments that are tidally disrupted, pushed all the way into 
the star (’’assimilated”), survive as gas giant planets at all separations, or survive as ’’warm giants” here defined as planets at separations 
less than 10 AU. 
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sible to notice a strong positive correlation with metallicity 
in the fraction of giants survived in the inner disc (’’warm 
giants”) and no such correlation in the outer disc. 


4 HOT PLANET MASS FUNCTION 

The planet mass function (PMF) is one of the most interest¬ 
ing observables of a planet population. Figure |3ji presented 
the mass function of planets formed in our simulations at all 
separations. The observed PMF are available only for close- 
in planets at this point in time as directly imaged planets 
are rare ('e.g.. lBowler et aDl2015l ~) and can only sample plan¬ 
ets more massive than ~ 1 Mj. We shall therefore focus our 
attention on the inner part of the mass-separation diagram. 

Since our model is still in the development stage, we 
shall not attempt a detailed comparison between the the¬ 
ory and the data by perfo rming a synthetic “observation” 
(see lMordasini et al.ll20095 l of our planets. Instead, we shall 
make a rather crude cut of the parameter space by carving 
out a “hot” planet sample which we define as the planets 
found within the inner 5 AU of the host star plus a 10% 
fraction of the ’’assimilated” planet group, that is, those 
planets that ended up at a = i?in. The latter group of plan¬ 
ets is included in our analysis for the reas ons explained in 
i|3.1lwh ich is also similar to the approach of iMordasini et ^ 
(l2009bl i. Specihcally, our inner boundary condition is too 
simple at this stage of the model development, yet observa¬ 
tions, especially transit surveys, are heavily biased towards 
these planets. We expect that some of the simulated planets 
at a = Rin will survive in the region between the star and 
a = 7?in. Our 10% choice of planets at a = Rin is ad hoc but 
this choice does not influence any of the conclusions of this 
paper significantly. We also consider below a planet sample 
that does not include the assimilated planets at all. 


4.1 Massive planets 

The top panel of hg.|T]shows the PMF of the hot planet sam¬ 
ple for planets with mass greater than 10 Mq . We show two 
versions of PMF. The histogram shown with the blue line 
shows all the planets in the sample. The shaded red color 
histogram shows the same sample but with a rudimentary 
RV-selection criterium of the stellar Dopier shift velocity 
exceeding ujet = 3 m/sec. To enable this selection, we gen¬ 
erate a random uniform distribution of cos i, where i is the 
inclination angle of the planet’s orbit (i = 0 corresponds to 
viewing the star-planet system along the rotation axis). The 
stellar Dopler’s velocity as seen by the observer is then 

V, fi^VKj^smi , (7) 


where vk = \JGMtja is the planet’s circular Keplerian ve¬ 
locity at its position, a, and M* = 1 Mq. We then prescribe 
the probability of detecting a planet with a given u* value 
as 


Pdet (^*) — 


V 


2 


Vf+V 


2 

det 


( 8 ) 


This is ad-hoc, but is asymptotically correct, so that all the 
planets with large v, 3> Vdet are detected, whereas none 
of the planets with v, <C Udet are detected. The smooth 
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Figure 4. Planet mass function (PMF) from simulations. The top 
panel is the “hot planet sample” defined in §4, with 10% of “hot 
Jupiter” planets included. The blue shows all planets whereas 
the red histogram shows planets affected by a velocity cut. The 
bottom panel shows planets in the same spatial disc region, with 
velocity cut imposed, but now excluding planets at a < 0.1 AU. 
The colours on the bottom panel separate metal rich and metal 
poor host populations. See E3] for detail. 


transition is more realistic than applying a sharp cutoff by 
requiring detection only for planets u* > Wdet- The mass 
shown in the histogram selected by w* > Wdet m/s is Mp sini 
rather than Mp, of course. 

The bottom panel of fig. 2] shows the PMF of the sim¬ 
ulated planets, with the velocity cut applied, and in which 
no planets inside the inner disc edge were included, so that 
a ^ 0.1 AU for all of these planets. We also separated the 
population on the metal rich (red) and metal poor (green) 
samples. 
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Atmospheres and Cores, a < 5 AU 



log Mass, 


Figure 5. Mass function of bound gas in the simulations (red 
histogram), compared with that of the cores (shaded green). The 
latter is scaled down by a factor of 3 for clarity. The bound gas is 
either self-bound (the high mass peak to the right) or is attached 
to numerous low mass cores (low mass part of the red histogram), 
with little in between. 


4-1.1 Tidal disruption desert 

One of the notable features of the simulated PMF, in both 
projections in £g. [4] is the deep depression in the number 
of planets with masses between ~ 20 Me and ~ 100 Me. 
An appropriate name for the feature is ’’Tidal Disruption 
desert” as it is due to the following. To the right of it, there 
are massive Mp ^ a few hundred Me gas-dominated plan¬ 
ets that did not experience a tidal disruption. To the left, on 
the other hand, are much smaller rocky core dominated plan¬ 
ets that are the most abundant remnants of the disruptions. 
The planets inside the desert, on the other hand, are the 
planets that went through a partial disruption. These plan¬ 
ets consist of a massive solid core dressed in a very dense 
gas envelope bound to the core as explained in section [2.2.3l 
The gas masses of these envelopes are comparable to or a 
few times larger than their core masses. We find that such 
planets are rare for the simulations presented in this paper. 

To expand on this point, fig. [5] compares the mass func¬ 
tion of all gaseous envelopes at Ri^ < a < 5 AU, shown 
in red colour, with that of the cores (shaded green his¬ 
togram). We see that gas envelopes (which may be quite 
metal rich, see © display a certain bimodality: they are 
either very low mass or very high mass. The low mass atmo¬ 
spheres are the majority; these are envelopes bound to low 
mass cores (Mcore a few M®). Higher mass cores attract 
more massive atmospheres, but then the number of massive 
(Mcore ^ 10 Mq) cores is small, and hence the number of 
more massive gas envelopes plunges as well. At gas mass of 
about 100 Mq « Mmin and above there is the domain of the 
undisrupted gas giants. 

While we are very confident that the Tidal Disruption 
desert is a robust feature of our model, its depth is probably 


not modelled reliably enough yet. The boun d atmosphere 
structure is dependent (iNavakshin et al.l |20l4) on the poorly 
know n opacity and the equation of state ( Stamenkovic et aH 
I2OI2I) of massive cores that are much hotter than present day 
cores of gas giant planets are believed to be. Additionally, 
our equation of state for the gas in the region very close 
to the core (see iHori fc Ikomal 20 111 , for how this may af¬ 
fect the critical core mass) is over-simplified. Formation of 
more massive atmospheres than found here would fill the 
desert somewhat, alt hough evaporation of atmospheres very 
close to the star (e.g., lOwen fc Wull20i:^ ). not modelled here, 
could empty it further. 

Note that Core Accretion model also predicts a 
“planet desert” at a sim ilar planet mass (llda fc Linll2004bl : 
iMordasini et al.ll2009bl) . for a seemingly different but actu¬ 
ally nearly identical physical reason. In CA, Mcore ~ 10 M@ 
is the critical core mass at which point a runaway accretion 
of gas onto the planet commences. Planets growing in the 
disc “travel” through the desert quickly to reach the higher 
mass (Mp rsj IMj) peak, so that the chance of a planet 
stranded inside the desert is relatively low. 

The lesson from this is that the structure of the planets 
may not always (or perhaps even rarely?) be used to con¬ 
strain the formation route of a planet (see also iHelled et al.l 
I2OI3I) . In the case of the desert in the PMF, atmospheres 
with masses ^ Mcore are unlikely, whether the planet is built 
bottom-up as in CA or top-down as in the TD. This is more 
of a fundamental property of matter than a consequence of 
a particular formation mechanism. It therefore appears that 
some other observational diagnostic must be used to distin¬ 
guish between CA and TD models. Metallicity correlations, 
dependences of the PMF on the mass of the host star, or 
timing of planet formation may hopefully break the degen¬ 
eracy in the future. 

The observed planet mass fu nction also s h ows a 
strong depression a bove ~ 20Mq iMavor et al.l 1 I 2 OIII) : 
iHoward et al.l (|2012h but it does not have the bump at ~ 1 
Jupiter masses that both our and CA models predict. We 
shall study this issue further in the future, but for now we 
note that a more efficient destruction of gas giant planets 
than our current calculations would reduce their number in 
fig- 51 s-nd also increase the number of planets in the ’’desert 
region”, making the models more compatible with observa¬ 
tions. 


4-1.2 Over-abundance of very massive planets 

Another feature of the theoretical PMF shown in the top 
panel of fig. 51 is that the number of very massive planets, 
Mp ^ 10 Mj « 3200 M 0 , is too large compared with the 
observed PMF. This region in the PMF is however sensitive 
to the initial fragment mass function shape, which we set to 
dNf/dMf oc 1/Mf (cf. fig. I3^). This would yield a flat PMF 
histogram. The distribution of initial fragment masses is not 
well constrained theoretically, and it may well be steeper 
than the one we used here. 

Furthermore, we also made the simplifying assump¬ 
tion (iv) in 112.11 that fragments do not accrete gas. 
It is possible that the most massive fragments do ac¬ 
crete gas rapidly (cf. simu l ations of INavakshin fc Cha 
I2OI3I : IStamatellos fc Herczed l2015l ) , so that they “run 
away” into the regime of low mass stellar companions 
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dStamatellos fc Whitworthll2008l ). which would remove them 
from our theoretical PMF. 

To expand on this, we note that the maximum 
available mass budget of gas to be accreted on a frag¬ 
ment is the mass of the disc, which is 0.15 Mq (see 
Table 1) for the simulations presented here. How¬ 
ever, this is a limitation of our set-up rather than a 
physical one. We here investigate only one fragment 
per disc, which makes most sense when the disc self¬ 
gravity is in its dying phase, and hence the disc mass 
is moderate. Real systems may however form gas 
fragments earlier on, in class 0 stage, when the pro¬ 
tostar is only a fraction of its central mass, and when 
the disc is still being fed from the protostar’s natal 
molecular envelope. If there is a “runaway” accre- 
tio n for more massive gas fragments as discussed 
in iNavakshin 8z Chal (l2013l l. then these fragments 
may grow to massive BDs and even low mass stellar 
companions. There is hence a potential connection 
of TD theory to low mass stellar companions, which 
may be another way to test it. This is however just a 
speculation in this point; gas accretion onto the frag¬ 
ments must be included in the calculations properly 
for reliable theoretical predictions. 

4.2 Core-dominated planets 

We now focus on core-dominated planets. We select the plan¬ 
ets less massive than Mp < 10^ ® M@, slice the resulting sub¬ 
sample into three spatial bins, and plot their PMF in figure 
[6l The histogram in blue shows the population of cores clos¬ 
est to the host star, in a < 0.5 AU region. The red histogram 
covers the region 0.5 < a < 5 AU, and the green histogram 
shows cores located further out. 


4-2.1 No pool of close-in very low mass cores 

Fig. [6] shows that the inner 0.5 AU region contains very few 
cores less massive than 1 M®. The absence of sub-Earths is 
due to the fact that tidal disruptions producing such low 
mass cores typically occur further out, at a few to 10 AU 
region. The low mass cores migrate slowly since the migra¬ 
tion rate is proportional to the core’s mass, cf. eq. [T] Most 
of them therefore do not make it inside the inner 0.5 AU 
region. 

This result will probably be weakened to some degree 
by more realistic calculations that would include previous 
generations of fragments during earlier disc phases. Cores 
locat ed in the inner few AU ma y get locked into resonances 
fe.g.. IPaardekooper et al1l2013tl with gas fragments migrat¬ 
ing from the outer disc and hence be pushed inward along 
with the fragment more rapidly than they could migrate 
by themselves. If the fragment is then tidally disrupted and 
the gas is consumed by the star, a second low-mass rocky 
core will be left behind. If observed at the present time, the 
system would have no obvious record of a past outer mas¬ 
sive fragment’s existence. Additionally, once the gas disc is 
removed, the N-body interactions of cores with each other 
may bring more low mass cores into the innermost region. 

Despite this, the rollover of the core mass function at 
low masses is a general result of our model for all separations 


(cf. fig. Ill), although the core mass at which the rollover oc¬ 
curs decreases with increasing separation. Therefore, while 
the exact core mass at the rollover is subject to further mod¬ 
elling, its existence is a robust result. 

We emphasise the difference from the Core Accretion 
way of building cores here. In CA, all of the solid mass is 
initially in tiny bodies (planetesimals). In the early stages of 
planetary growth, most of the mass is in the low mass cores 
and “embryos” with mass well below 1 M@. The PMF from 
CA calculations thus shows a strong peak at low masses, 
M 1 Mffi, see fig. 3 iniMordasini et al.l ll2009bl l and fig. 2 in 
iMordasini et al.l (l2012ll . There is no physical reason to have 
such a divergence towards smaller masses in TD theory for 
planet formation since the cores grow by accretion of grains 
and not by collisions of numerous solid “embryos”. 

We note that this pr ediction of the mod el is consistent 
with Kepler observations llHoward et al.l[2bl2l') that find that 
the frequency of planet occurrence (corrected for observa¬ 
tional biases) does not diverge and instead drops towards 
the smallest radius planets {Rp ~ R©). 

4-2.2 A smooth PMF with no local minima at a few M© 

Another interesting distinction between our calculations and 
that of CA theory is that some of the latter predict that 
the cores can be composed of two separate populations 
- the inner smaller rocky cores and the outer more mas¬ 
sive icy cores. Some of CA models (e.g., see figs. 8 and 
13 iniMordasini et al.l (l2009al ) and fig. 3 in iMordasini et al.l 
(l2009bl u show a significant dip in the PMF in-between these 
two populations, at Mcoie ~ a few M©. The fact that the 
low mass cores are rocky and the high mass ones are icy is 
best s een in figures (especially 5, 7 and 11) in lAlibert et al.l 

1I2OI3II . 

As will be seen later, our cores, including the most 
massive ones, are rock-dominated, and therefore there is no 
physical reason for two separate populations to exist in the 
PMF of the cores. It is hence a smooth function of Mcore 

(fig-®- 

4 . 2.3 A rollover above ~ 10 M© 

The population of cores shown in fig. [6] in the inner 0.5 AU 
is dominated by super Earths of mass 5M©. There are 
very few cores more massive that 10 M©. As noted in 1)2.2.21 
in TD cores grow by grain sedimentation, and this process is 
limited by the rate at which the grains can sediment without 
fragmenting in high speed collisions, and also by the convec¬ 
tive mixing of the grains which tends to bring the grains from 
the centre bac k into the outer regions o f the fragment (cf. pa¬ 
per I and also lHelled fc Bodenheirn^ ll2010li '). Massive cores 
are very luminous in our model. This amplifies convection in 
the inner parts of the fragment, forming a kind of a feedback 
loop that limits the rate at which the cores grow. 

Therefore, the absence of more massive cores in our 
models appears to be due to an insufficient time during 
which the conditions within the fragments are conducive for 
core growth. For a rapid core growth, the fragment must 
be (a) dense enough to allow an efficient grain growth and 
sedimentation, (b) but not too hot for the grains to vapor¬ 
ise and for the fragment to collapse into the second core 
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Simulated cores 



Figure 6. Planet mass function (PMF) of cores from simula¬ 
tions in three different spatial regions as explained in the legend. 
The vertical scale on the right shows the cumulative probability 
distribution for the same curves. See EJ for detail. 


Simulated Planets, a < 0.5 AU 



when the core will be buried under all the weight of the gas. 
The core growth time window is also limited by the lifetime 
of the fragment from its birth to the moment it is tidally 
disrupted. 

The precise value of the upper rollover mass in the core- 
dominated part of the PMF does depend on parameters of 
the models and the atmosphere presence for the most mas¬ 
sive cores. To explore this issue, we slice the sample of plan¬ 
ets for a < 0.5 AU into two sub-population for which the 
grain breaking velocity, Ubr, is between 5 < Ubr < 7 m s~^ 
and 12 < Ubr < 15 m s“^, and plot the resulting mass func¬ 
tion in fig. [T] The higher breaking velocity sub-population 
naturally has more massive cores because the grains can sed¬ 
iment at higher velocity. This shows that the exact shape of 
the mass function near Mp ~ 10 — 30 M® is not well con¬ 
strained at the moment. Additionally, it is possible that ac¬ 
cretion of gas or pebbles after the tidal disruption of the 
fragment would increase the mass of the planets shown in 
fig. [S] somewhat. 

The red histogram in fig. [^shows that low mass planets 
(mainly bare cores) become more dominant at larger sepa¬ 
rations. This is due to the fact that high mass cores migrate 
inward much more than do low mass cores. The green his¬ 
togram indicates that there are relatively few cores beyond 
5 AU. This result is however sensitive to the migration pre¬ 
scription and the speed with which the most massive cores 
are assembled; fig. 4 in paper II shows a much larger popu¬ 
lation of cores formed at tens of AU from the host star. 


5 PLANET METALLICITY PREFERENCES 

We shall now analyse how the host star metallicity influ¬ 
ences the likelihood of planet formation of different types. 
To reiterate the setup described in 112.2.41 our metallic¬ 
ity distribution is a Gaussian (equation [S]) with the mean 
[Fe/H] = [M/H] equal to zero and the dispersion around 
that a = 0.22, motivated by the CORALIE planet sample. 
The nature of the resulting dependences turn out to vary 
from a region to a region in the planet’s mass-separation 
plane; thus we shall study these regions separately. 


Figure 7. Mass function of low mass planets in the inner 0.5 
AU for low and high grain breaking velocities, as specified in the 
legend (ubr is in units of m/s). Unsurprisingly, more massive cores 
are assembled within the fragments if grains can sediment without 
fragmenting at higher velocities. 


5.1 Moderately massive giants vs super Earths in 
the inner 5 AU 

We now analyse the metallicity correlations in the hot planet 
sample (see 114.11 the planets found within the inner 5 AU of 
the host star plus a 10% fraction of the ’’assimilated” planet 
group, that is, those planets that ended ended up at a = 
Pin). From this sample we then select super Earths defined 
as planets in the mass range 2M@ < Mp < 15 Mq, and 
moderately massive giants as gas giants with mass between 
at least 100 M0 « 0.3 Mj and 5 Mj. More massive giants are 
rare in the observed samples. There is also a physical reason 
to look at the metallicity correlations of the more massive 
giants separately, as we will see in 115.31 

Fig. |8] shows how these two groups of planets are dis¬ 
tributed over the host star metallicities. The blue line his¬ 
togram shows the giants, whereas the red histogram shows 
the super Earths. The curves of the same colour show the re¬ 
spective cumulative distributions with the scale on the right 
vertical axis. If our theory were insensitive to the metallicity 
of the host star, then we would expect to see a Gaussian- 
like distributions centred on zero in fig. [5] However, it is 
evident that gas giants appear preferentially around metal- 
rich hosts, whereas super Earths are spread about the mean 
metallicity roughly symmetrically. Also note that the fre¬ 
quency of finding a gas giant per star of a given metallicity 
keeps increasing with [M/H] up to the highest value of 0.5 
shown in the figure (see fig. [TU] below). The reason why the 
gas giant histogram goes down in fig. [8] is that there are very 
few stars with meta llicities larger than [M/H]=0.3. 

As explained in iNavakshirJ ll2015bli and papers I and H, 
the positive metallicity correlation for giants is due to peb¬ 
ble accretion rate onto the fragment being larger at higher 
metallicities, which leads to the fragment contracting faster. 














































Simulated Planets, hot region 
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Figure 8. The fractional distribution of gas giant planets (blue 
filled-in histogram) and super Earths (red histogram) over metal- 
licities of the host stars. Only planets in the inner 5 AU are se¬ 
lected. The curves of the respective colours show the cumulative 
distributions of same planet groups. Note strong positive metallic- 
ity correlation for gas giants and absence of such for super Earths. 
See text in ani for more detail. 


Faster contraction implies fewer tidal disruptions, and hence 
more gas giants surviving at higher metallicities. The phys¬ 
ical reason for pebble accretion being so influential in con¬ 
trolling the rate of the fragment’s contraction is that pebbles 
bring in additional mass into the fragment (but not kinetic 
energy since they sediment onto the cloud gently or else they 
break up in high speed collisions). Gas clumps dominated by 
molecular hydrogen turn out to be very sensitive to addition 
of mass in this way and collapse when the mass of the extra 
metals reaches ~ 5% to ~ 20% of the initial fragment mass. 
This picture also has implications for metal overabundance 
in giant planets as a function of their mass as we shall see 
in 1|S] The more massive the fragment is, the hotter it is 
at birth, the less pebbles is required to bring it to the H2 
dissociation instability point (when the central fragment’s 
temperature is « 2000 K). 

The metallicity correlation for super Earths turns out to 
be a much subtler matter. Fig.[8]shows that high metallicity 
environments are no more preferable for formation of cores 
in the mass range 2 Mq < Mp < 15 M0 than low metallic¬ 
ity discs. In paper II, two explanations for this result, both 
working in the same direction, were proposed: (a) saturation 
in the mass of the most massive cores within the fragment 
at high metallicities, and (b) disruption of high metallicity 
gas fragments being too rare. The latter effect is important 
because super Earths are the remnants of these disruptions: 
no disruptions, no super-Earths. 

We now investigate the relative importance of these two 
effects. The top panel of Fig. [9] shows the mass function of 
cores found in our simulation in the inner 5 AU for four 
metallicity bins as specified in the legend. This panel shows 
that most of the super Earths form at metallicity not too 
dissimilar from the Solar metallicity (the blue and the green 
histograms). The curves are normalised by the total integral 
giving the total number of cores, of course. The bottom panel 
of the hgure shows the distribution of frequency of planet 
masses for the curves in the top panel; the histograms on the 
bottom are normalised on unity for each metallicity bin. It 
is clear from the bottom panel of hg. |9]that the mean mass 
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Simulated cores, a < 5 AU 
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Figure 9. Number (top panel) and fractional (lower panel) mass 
distributions for cores found in the inner 5 AU for four different 
metallicity groups. 


of the core increases with metallicity rapidly while [M/H] < 
0.25. However, the mean core mass increases less between the 
two metal rich bins. There is thus indeed a certain saturation 
in the cores mass growth at high metallicities, although the 
mean core mass still increases with [M/H] . Thus effect (a) 
is not actually that strong. 

The main driver of the poor correlation between super 
Earths and host’s metallicity in our models turns out to be 
(b), e.g., that the number of tidal disruptions is just too low 
at high [M/H], yielding too few cores (although on average 
they are more massive than the cores at lower metallicities). 
Figure [To] shows the frequency of the fragment becoming a 
given planet type as shown in the legend versus metallicity 
of the host star. In particular, the black connected diamonds 
show the frequency of obtaining a super Earth in the inner 
disc. The red crosses show the frequency of the fragment 
migrating through all of the disc to r = i?in without be¬ 
ing disrupted. Notice the anti-correlation between such un¬ 
disrupted fragments and the number of super Earths. In par¬ 
ticular, at the highest metallicities, as much as ~ 70 — 80% 
of our initial fragments manage to collapse and migrate all 
the way in, avoiding tidal disruption. This of course must 
mean that fewer core-dominated planets are made in this 
environment because there are fewer tidal disruptions. 

Summarising these findings, we conclude that Super- 
Earths do not correlate with metal abundance of host star 
in our models because gas giant planets do. The presence of a 
gas giant planet indicates that its predecessor, a pre-collapse 
molecular fragment, was successful in avoiding the tidal dis¬ 
ruption. Since gas giants are most frequent at high metal¬ 
licities, gas fragment disruptions must be much less ubiqui¬ 
tous at high [M/H], and hence there should be fewer Super- 
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Cores and low mass giants vs [M/H] 



Simulated Planets 



Figure 11. Same as fig. [8] but now comparing the metallicity 
distribution of gas giant planets in the inner 5 AU (hot giants) 
with those in the "cold” region, a > 10 AU (red histogram). Note 
that cold giants show no metallicity preference. 


Figure 10. Frequency of a fragment becoming a planet of a given 
type, as shown in the legend, as a function of the host star metal¬ 
licity. Note that at the highest metallicity bins, most fragments 
are able to migrate in to a = R[n and hence only a few super 
Earths are made. Also note that it is the blue square line that 
should be compared with Fischer Sz Valenti (2005), since this line 
shows the frequency of finding a gas giant per star of a given 
metallicity, rather than that in the stellar population as a whole, 
which is presented in fig. 8. 


Earths. On the other hand, we note that the mean mass of 
the cores does increase with metallicity. Thus, Super-Earths 
made by TD model may be said to correlate in mass (some¬ 
what weakly at high core masses) but not in numbers with 
the metallicity of the host star. 

The histogram distribution for our simulated planets 
shown in fig. [8] looks quantitatively similar to fig. 2 in 
iBuchhave et aU (|2012l )' even though we did not make any 
attempt to arrive at such an agreement specifically. 


5.2 Cold giants are metal insensitive 

We now compare the metallicity correlation for all of the 
gas giant planets Mp > 100 M 0 in the hot sample of i|4.1l 
with the same mass range planets located further out, at 
7? > 10 AU, which we name “cold giants”. 

Figure ITT] compares the metallicity distributions of the 
hot and cold giants. As was the case with Super-Earths dis¬ 
cussed in il5.ll the cold giants show little metallicity prefer¬ 
ence. The mean of the hot giant distribution is at [M/H] = 
0 .12, clearly showing a significant positive metallicity pref¬ 
erence, whereas for the cold giants we find mean[M/H] = 
-0.01. 

The interpretation of this result is straightforward. 
Most of the fragments that stayed behind 10 AU did not 
come close enough to the parent star to experience tidal 
forces sufficiently strong to rip the fragments apart. There¬ 
fore, these planets survive irrespectively of how much peb¬ 
bles was accreted. They are therefore insensitive to the metal 
content of the parent discs. 


5.3 Most massive giants are metal insensitive 

Finally, we come to the most massive planets studied in 
this paper. A visual inspection of the top left hand corner of 
Figure[2]shows that metal poor massive giants (Mp ~ 10 Mj) 
are about as abundant as metal rich planets of same mass. 
This appears true for planets at the inner disc edge and also 
in the inner few AU, so it may be general for the whole hot 
sample of the most massive planets. 

Figure [ 12 ] shows that this is indeed the case. Similar to 
fig. HOI this figure shows how frequency of an initial frag¬ 
ment becoming a given type of planet changes with the host 
star metallicity. Three groups of gas giants is considered in 
fig. [T 2 I The fraction of survived moderately massive giants 
(the blue diamonds curve) is a rapidly increasing function 
(roughly a power law oc where Z = re¬ 

flecting the strong positive correlation for the hot moder¬ 
ately massive giants discussed in il5.ll 

The black squares curve in fig. [T^ shows the survival 
fraction for more massive giants, Mp > 5Mj. We see that 
the probability of these massive fragments surviving the disc 
migration phase inside the hot region is actually flat with 
metallicity. Further to that, the red crosses show the fraction 
of the initial fragments in the massive sample that migrated 
all the way to the innermost disc radius. The curve has a 
depression at metallicities just below the Solar value, which 
is to say that fragments around stars of nearly Solar compo¬ 
sition are the most likely to be tidally disrupted before they 
arrive at a = Rin- 

For the most massive fragments, therefore, very high or 
very low metallicity environments are actually preferable (al¬ 
though not by much, e.g., less than by a factor of 2). This bi¬ 
modality of metallicity preferences of hi gh mass giants form - 
ing in the context of TD was hinted on in iNavakshinI (l2015al f , 
where it was shown that gas giant planets contract and col¬ 
lapse in one of two ways. The most well know n way, studied 
in lite r ature on giant plan e t evol ution since iBodenheimerj 
||1974| : iBodenheimer et ^ lll980lj . is by radiative cooling, 
when the planet radiates its excess thermal energy away be¬ 
fore collapsing. The second way is by metal loading via peb¬ 
ble accretion, in which case it is the increasing weight that 
leads the planet to collapse. The radiative cooling time of the 
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Massive giants vs [M/H] 



[M/H] 

Figure 12. Similar to fig. 1101 frequency of the initial fragments 
becoming a planet in the three mass-separation ranges defined 
in the legend. Note that low to moderate mass giants (blue line) 
correlate positively with metallicity, whereas high mass giants do 
not. 


fragments is a str ongly decreasing f unction of the fragment’s 
mass (cf fig. 1 in INavakshinll^lSal ) . In section 7 of that pa¬ 
per, it was suggested that even at large grain opacities, that 
is not modified by grain growth, high mass fragments are 
capable of contracting faster by radiative cooling than by 
pebble accretion. 

The radiative cooling rate of the fragment decreases 
with grain opacity, which is directly proportional in our 
model to the total mass of metals in the planet minus the 
mass of the core. Therefore, at higher metallicities the radia¬ 
tive cooling channel for collapse of massive fragments is no 
longer efficient. However, pebble accretion does help at the 
highest metallicities, explaining why the curve shown with 
the red crosses recovers to almost unity at high [M/H]. 

Summarising this discussion, the metallicity correlation 
for the highest mass giant planets is essentially a sum of two 
correlations: (i) an anti-corr elation for the radiative l y cool - 
ing fragments predicted by iHelled fc Bodenheimeil (1201 ll ) 
and (ii) the posit ive correlation due to pebble accretion 
llNavakshirJl2015^ ). The reason why the anti-correlation be¬ 
comes important only for high mass Mp ^ 5 Mj planets and 
not lower mass planets is that for lower mass planets the 
radiative cooling channel is not important at all compared 
with the pebble accretion, even at very low metallicities. 

One caveat must be emphasised here. We do not include 
accretion of gas onto pre- or post-collapse giant planets (see 
point (iv) in i)2.1ll . It is conceivable that some of the ob¬ 
served massive Mp ^ 5 Mj planets could have started from 
gas fragments of lower mass (even in the context of TD). 
If this is the case then those planets would correlate with 
metallicity positively since their progenitor fragments did. 
This would dilute the effects we discussed here, but we are 
not yet able to quantify by how much. 


Metallicity correlations in the Mass vs Separation plane 



Figure 13. The planet mass versus host separation plane, with 
l/40th of the simulated planets shown, shaded by the type of 
metallicity correlation in the respective region. “Number positive” 
and “Mass positive” are correlations in which the frequency of 
planets or the only the mass of the planet correlates with the 
star’s metallicity, respectively. See text in 3131 for explanation. 
The symbol colour meanings are same as in fig. [2] 

5.4 Metallicity preferences, simulation vs data 

In the previous sections we presented theoretical predictions 
on how planets correlate with the host star’s metallicity in a 
number of regions in the planet mass - host star separation 
plane. For convenience, these are summarised graphically 
in fig. [T3] in which we show a small subset of the simu¬ 
lated planets (l/40th fraction of the total, uniformly and 
randomly selected for the whole plane) overlay-ed by three 
shaded regions. These are: 

(i) In the top green shaded region, planets are insensitive 
to the host star’s metallicity. For cold giants ( >15.211 this is 
because these planets do not experience strong tidal forces 
at their far away locations. They hence eventually collapse 
independently of the presence of pebbles in the host disc. For 
hot very massive planets, a second collapse route - by ra¬ 
diative cooling ~ becomes available, which “saves” the frag¬ 
ments from tidal disruptions at low metallicities ( >15.311 . 

(ii) In the red region, pebble accretion is most effective 
at preventing tidal disruptions at high metallicities, so there 
is a strong positive metallicity correlation ( >15.111 . Although 
we did not analyse this specifically, there is also a positive 
metallicity correlation for partially disrupted planets, e.g., 
with masses approximately between that of Neptune and 
Saturn, because these planets have very massive cores, and 
massive cores are assembled preferentially at high metallic¬ 
ity. 

(iii) In the black shaded region, the metallicity correlation 
is complicated. If measured within a certain mass window, 
there is no clear correlation for the number of cores ver¬ 
sus [M/H]. For example, cores with mass Mcore ~ SMq are 
relatively rare at low metallicities (since in this case most 
cores are less massive than that, see fig. O, then abundant 
at ~ Solar metallicity, and are rare again at high metallici¬ 
ties, since most cores then become more massive than 3 Mq 
. However, there is an average core mass correlation with 
metallicities: the higher the metallicity, the more massive 
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Simulated Planets, a < 5 AU 
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Observed planets, a < 5 AU 
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Figure 14. Mass function of simulated (top) and observed (mid¬ 
dle) planets, sliced into the metal rich and metal poor samples. 
The bottom panel shows the fraction of planets that are metal 
poor for the three selected groups of planets, for both observa¬ 
tion and simulation. No attempt to include errors or observational 
biases in this figure was made. 


an average core. We called this type of correlation ’’number 
neutral but mass positive” in gsu 

It is tempting to compare these theoretical predictions 
to the observations. Observations of planets in the cold re¬ 
gion are too incomplete for us to attempt a comparison, 
hence we limit our attention to the inner 5 AU region of 
the simulated planets. As we wish to compare this to ob¬ 
servations, we included only planets with the stellar reflex 
velocity u* > 2 m/s. This preferentially removes the planets 
less massive than ~ tens of M 0 from the simulated sample. 

We have no access to a uniformly selected observed 
planet sample for the broad mass range to investigate metal- 
licity correlations of the observed planets. Therefore, we only 
attempt a rather simple comparison here. We assume that 
metal-rich and metal-poor planets of same mass are affected 
roughly similarly by detec tion biases, and use th e planets in 
the exoplanets.eu sample (ISchneider et al.| 1201 il l. We select 
planets with measured masses and host-star separations less 
than 5 AU, and with the host star masses between 0.6 Mq 
and 1.5 Mq. 

We slice the populations of the simulated and observed 
planets into the metal poor ([M/H] < 0) and the metal 
rich ([M/H] > 0) samples. In fig. 1141 the top panel shows 


the PMFs for the simulated planets inside the inner 5 AU. 
The blue, the red and the hashed green histograms show 
the PMF for the total, the metal rich and the metal poor 
samples, respectively. The middle panel in he 1141 shows the 
observed sample of planets (cf. 6g. [4]), also sliced by the 
metallicity of the host. No attempt to include errors or ob¬ 
servational biases in this panel was made. For this reason, 
the observed PMF shown in the bottom panel massively 
under-estimates the number of low mass planets which are 
much harder to detect than high mass ones. However, the 
ratio of the metal poor to metal rich systems at a given 
planet mass is probably much more robust, and this is what 
we indent to compare with the simulations. 

The bottom panel of 6g. [14] shows the fraction of stars 
that are metal poor for three mass groups of planets. Qual¬ 
itatively, similar trends are seen in the simulations and in 
the observations: the Super-Earths are (number) insensitive 
to the host star metallicity, medium mass giants are posi¬ 
tively correlated with the host star metallicity, and the most 
massive gas giants are somewhat weakly correlated with the 
metallicity. The highest mass bin in this picture is particu¬ 
larly interesting since, it appears to us, in CA the highest 
mass giants should still correlate with the metallicity. These 
planets are formed in CA by accreting more gas onto less 
massive giants, and those do correlate with z. There is thus 
a potential to distinguish CA and TD by the metallicity 
correlations for the highest mass gas giants. 

As already mentioned, the observational sample for ex¬ 
oplanets used here is unfortunately not u niformly selected 
and m ay have various selection biases. lAdibekvan et all 
ll2013ll presented a study of metallicity correlations for plan¬ 
ets selected in a much more homogeneous way from the 
SWEET-Cat database. Among a number of correlations, 
they find that most massive giants in their sample (more 
massive than 4Mj) correlate with metallicity weaker than 
less massive giants, although the statistical signihcance of 
this result is not very robust due to a small number of plan¬ 
ets in the high mass bin. 


6 ATMOSPHERES ON TOP OF CORES 

As described in 1)2.2.31 in this paper we also calculate the 
mass of gas gravitationally bound to the cores inside the 
pre-collapse fragments. When the fragment is disrupted this 
atmosphere is assumed to survive, remaining bound to the 
core. This is a reasonable assumption, since the density 
of these atmospheres are typically higher than that of the 
main body of the frag ment by orders of magnitude (e.g., cf. 
iNavakshin et al.1l20l4l l. The atmospheres could loose more 
mass after the fragment disruption by photo-evaporation 
if the remnant mi grates sufficiently close to the star (e.g., 
lOwen fc WrJ[2oiT : this post-disruption evaporation is not 
included in our calculations at this stage. 

The top panel of figure [15] shows the ratio of the atmo¬ 
sphere mass to that of the core for the simulated planets in 
the inner 5 AU. Colours and types of the symbols indicate 
the host star metallicity, as before. We see that metallicity 
does not play a big role in establishing the mass of the atmo¬ 
sphere at a given Mcore. The metallicity is more important 
in determining the core mass itself, as we saw in i)5.1l One 
notices that metal rich systems (blue and red) appear more 
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often at Mcore ~ SMq part of the figure, whereas the metal 
poor systems typically have smaller mass cores. 

Notice that there is a large spread in the relative atmo¬ 
sphere mass at a given Mcore- This spread is not a function 
of metallicity. The physical reason for the spread is that 
cores of a given mass can be born inside fragments of differ¬ 
ent masses and/or different evolution histories. These frag¬ 
ments have different central gas pressure and densities which 
explains why same mass cores may have different atmo- 
sph ere masses as measure d at the fragment disruption (see 
also lNavakshin et al.ll201^ i. In addition, we parameterise the 
core’s luminosity by introducing the Kelvin-Helmholtz time 
for the core, tkh, which is a Monte Carlo parameter varied 
between 10® and 10^ years (see Table 1 and 1)2.2.31) . Due to 
this, the luminosity of the cores vary, and that also imprints 
onto the scatter of the atmosphere’s mass at a given Mcore- 

With all the scatter, it is nonetheless clear that cores 
less massive than ~ 5 Mq do not possess atmospheres more 
massive than a few percent of the core, whereas cores more 
massive than that may have atmospheres of mass compa¬ 
rable to that of the core. A few of these formally had at¬ 
mospheres more massive than the core, which we capped 
at Matm = Mcore exactly as explained in 1)2.2.31 Physi¬ 
cally in these cases we expect that a tidal disruption of the 
fragment will result in a partial disruption of the fragment 
only as there can be more mass in the collapsed atmosphere 
around the core than the core mass itself. At the same time, 
we note that in fig. [15] there are also quite massive core, 
Mcore (i/ 10 M 0 that have very small, 1 %^ atmospheres. 

The bottom panel of fig. [15] shows the metallicity of 
the atmospheres of the cores, -Zatm- All of these values are 
strongly super-Solar {Zq = 0.015). This is to be expected as 
gas fragments accrete pebbles in our model, so their abun¬ 
dances are super-Solar (cf. 3H] on this for more detail). In 
addition to that, grains do sediment to the centre of the 
fragment and so the metal abundance in the atmosphere 
is usually even higher than the mean fragment metallicity 
(e.g., figure 9 in paper I shows that the metallicity is the 
highest in the centre of a planet). 


7 ROCK DOMINATED COMPOSITION OF 

CORES 

Figure [in] shows the ratio of the volatile (water and CHON) 
mass in the core to the total mass of the core, again for 
the four different metallicity groups. We see that with a few 
exceptions at the low core mass end, the cores are always 
dominated by the rocks rather than “ices”. This is especially 
true at the high core mass end, where CHON contribute a 
few % of the mass only. 

This result is not particularly surprising in the light 
of previous studies of grain sedimenta t ion in side the pre¬ 
collapse fragments. iHelled fc SchubertI (l2008l 'l showed that 
water ice grains are unlikely to sediment into the cores of 
pre-collapse gas fragments of mass higher than about 1 Mj 
because the fragmen ts are too hot for th e icy grains to exist 
even at early times. iHelled et al.l l|2008l ') then showed that 
organics (modelled as CHON) grains have somewhat better 
chances of sedimenting into the cores since their vaporisation 
temperatures are 350 to 400 K, depending on gas pressure, 
whereas water ice vaporises above T ~ 150 — 200 K. Rocks 
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Figure 15. Top: the ratio of the mass of the bound atmosphere 
to that of the core as a function of core mass after the fragment 
disruption. Bottom: the metallicity of the atmospheres shown in 
the top panel. The legend explains the meaning of the symbols. 


and Fe (counted together as one species in our paper) are 
thus the main compon ents of which the cor es are made in 
the TD model (see also lForgan fc Ricdl2013l ’). 

There remain quantitative uncertainties in the compo¬ 
sition of the cores. We here chose a rather low initial tem¬ 
perature for the fragments, so that water ice can sediment 
initially in gas fragments less massive than a few Jupiter 
masses. This is why there are some volatile grains in the rel¬ 
atively low mass cores in fig. 1161 However, high mass cores 
require a long, a few 10“* years, assembly time. The frag¬ 
ments containing those heat up to temperatures of hundreds 
of K during this time, so that neither CHON or water ice 
grains can sediment. 

This prediction of the model is important as it is very 
different from CA model in which massive cores have the 
best chance of growing beyond the ice line where water ice 
can sediment, so that they are dominated by ices. 


8 METAL OVERABUNDANCE OF GAS 
GIANTS COMPOSITION 

Gas giant planets, and also their less massive cousins, 
are over-abundant in elements heavier than H/He com- 
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Figure 16. The ratio of the core mass in the volatiles (water ice 
and CHON) to the total mass of the core as a function of the core 
mass. Note that high mass cores are strongly dominated by rocks 
in our model, in contrast to CA. 


Figure 17. Ratio of metal abundance of the giant planets to that 
of their host star. The symbols with the error bars are models by 
Miller Sz Fortney 2011 based on observations of close-in but not 
strongly irradiated gas giant planets. 


pared with their host sta rs. This is well kn own for the 
Solar System planets (e.g.. [Ciiillot et ^l2004l and is now 
also robustly establishe d for exoplanets at close separations 
dMiller &: Fortnev 20 111) . This important fact is usually used 
(e.g.. iHelled et J 120131) to argue that the CA model is a 
clear favourite to explain giant planets formation because 
planets made by CA contain heavy cores which makes the 
planets more compact (as required by the observations) than 
they would be if they had their host star metallicities and 
did not contain massive cores. In contrast, GI planets may 
be enrich ed in metals by accretio n of grains or planetesimals 
early on (|Bolev fc Duriser]l201(Th . but until recently it was 
thought that this enrichment is not fundamental to the sur¬ 
vival of GI planets. It would thus seem that GI planets may 
have a range of metallicities - fro m sub-host values to many 
times the metallicity of the host llBo ley et aO 2^11^ . 

However, as argued in iNavakshinI ( 2015al l. pebble ac¬ 
cretion is the process that makes relatively low mass (a few 
Mj) GI fragments survival possible. Radiative cooling is too 
inefficient, and it is pebble accretion that allows these frag¬ 
ments to contract into the much denser “hot start” GI planet 
configurations so that they avoid tidally disruptions in the 
inner disc. Therefore, we may expect that the metallicity of 
TD giant planets is elevated compared to their host stars. 

Figure [17] shows the ratio of the metallicity of the simu¬ 
lated planets in the inner 5 AU, Zpi, to the metallicity of the 
host star as a function of the planet’s mass. As before, the 
population is broken into four bins by the host star metal¬ 
licity, [M/H], as indicated in the legend. I n addition to that, 
we us ed the “mixed model” in Table 1 of iMiller fc FortnevI 
ll201l]i to show their best estimates for such a ratio for their 
sample of observed giant planets. These models are shown 
with the blue symbols with the error bars in fig. 1171 

Fig. indicates that TD giant planets are indeed 
over-abundant in metals compared to their host stars. In 
addition, there is an anti-correlation between the “over- 
metallicity” ratio, Zpi/Zstar, and the mass of the planet. The 
origin of this anti-correlation is slightly different for plan¬ 


ets that did not experience tidal disruptions, Mp ^ 0.33 Mj 
here, and those that did (the less massive ones). 

For the partially disrupted planets, the correlation is 
physically similar to that of GA gas giant planets. The metal 
content of these planets is actually dominated by the core’s 
mass, which is Mcore ~ 10 Mq. The atmosphere, while also 
over-abundant in metals as we saw in fig. 1151 contains a 
minority of the metals mass. Since M^ore ~ const, and is 
independent of the total planet mass, Zpi ~ is this 

regime. 

For heavier planets that did not go through a tidal dis¬ 
ruption, the correlation occurs bec ause (a) higher mass plan¬ 
ets are hotter to begin with (e.g.. iHelled fc Schubertll2008l : 
lNavakshI^l201lh and (b) collapse faster as radiative coolin g 
contributes more to their contraction (INavakshinI l2015a] l. 
Therefore, more massive gas fragments require a smaller 
fractional mass increase in pebbles to reach the central tem¬ 
perature of ~ 2000 K and collapse. 

Note that we do not take into account gas accretion onto 
the planets in this paper in either pre- or post-collapse evo¬ 
lution. If this assumption was relaxed, post-collapse planets 
could gain mass by accretion of gas with metallicity close 
to that of the host star. This would therefore dilute their 
metallicity at a given total planet mass, but, importantly, 
would not change the downward trend that we see in fig. 1171 

Our simu lation results are qualita tively consistent with 
the models of iMiller fc FortnevI 320111 '). although we did not 
fine tune any of the parameters to specifically address this is¬ 
sue. We finish this section reiterating what we said in i)4.1.1l 
the present day bulk structure of the planets may not be 
uniquely indicative of the route that these planets formed. 
We believe TD is as promising as CA in explaining the bulk 
composition of Solar and extra solar planets. 
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Planet vs fragment mass 



1 10 
Fragment mass [Mj„p] 



Fragment mass [Mj^p] 


Figure 18. Top: Final mass of the planet versus the initial mass 
of the fragment. The meanings of the symbols is the same as in fig. 
1151 Bottom: Mean mass of the planets more massive than 5 M0 
versus initial fragment mass for three regions in the host-planet 
separations. 


9 WHICH FRAGMENTS MAKE WHICH 

PLANETS? 

One may wonder how the value of the initial fragment mass 
maps into the kind of planets that emerge in the end. To 
investigate that, the top panel of figure [18] shows the final 
mass of the planet versus the initial fragment mass, Mf, for 
a uniformly randomly chosen (7% of the total) subsample 
of all the planets in the simulations. The colours reflect the 
metallicity of the system, as before. The nearly straight line 
in the upper part of the panel are the fragments that avoided 
tidal disruption and collapsed into giant gas planets. Mass 
of these planets is close to the initial fragment mass, save 
for the addition of metals by pebble accretion; hence the 
(almost) one to one correspondence between the final and 
the initial planet masses. 

The planets below that line are planets that went 
through a total or a partial gas envelope disruption. These 
are mainly rocky core planets with only a few systems that 
are dominated by gas, as we saw earlier. 

This figure demonstrates that the mass of the core made 
within the fragment is not a unique function of the frag¬ 


ment’s mass. In fact, the most massive cores are found in 
the least massive f ragments, those bel o w Mf ^ 2Mj. This 
result is not new. I Helled et al.l (|2008l i: iHelled fc SchubertI 
(l2008li showed that the more massive the fragment is, the 
less massive is the core made inside it during the pre¬ 
collapse phase. This is because the central temperature of 
pre-collapse fragment is believed to increase with its mass 
even at the initial stage due to adiab atic compression being 
stronger for more massiv e fragments (iMasunaga fc Inutsukal 
I 2 OOOI : NavakshinI I 2 OIOIJ 1 . Massive fragments also cool ra- 
diatively more rapidly, cutting the available time for the 
core’s growth. The exact value of the initial temperature, 
To, is however difficult to calculate from first principles, es¬ 
pecially in the case of the disc. In this paper we assumed 
that To oc My^. A stronger dependence would make even 
smaller cores at high Mj than the top panel of figure 1181 
shows. 

The bottom panel of fig. [T5] presents the mean mass of 
the planet as a function of the initial fragment mass. This is 
calculated neglecting the least massive planets, Mp < 5 M 0 , 
assuming that these are unlikely to be detected due to selec¬ 
tion bias. Further, three spatial regions are picked for this 
panel: the planets on the very edge of the disc (red crosses), 
planet found anywhere in the disc (blue squares) and planets 
found between the inner edge and 5 AU (black diamonds). 

In general, we see that the mean mass of the resulting 
planet correlates positively with the mass of the initial frag¬ 
ment. This result is of no clear importance if the spectrum of 
initial fragment masses is same for all the systems studied - 
as assumed in this work. This may not necessarily be so, and 
further work on the gravitational fragmentation of discs is 
needed to study this issue. One rather likely outcome, how¬ 
ever, is that the mass of the fragment would increase with 
the mass of the host star. If this is the case then the bot¬ 
tom panel of fig. [18] predicts that the mean mass of a planet 
would increase with the mass of the host star. 

The most massive cores, Mcore ~ 10 M 0 are made in 
relatively low mass fragments, M ^2 Mj, as emphasised be¬ 
fore. We can hence expect a similar result for massive atmo¬ 
spheres bound tho the core. Fig. 1191 indeed shows that cores 
with massive atmospheres are only made in low mass frag¬ 
ments. This result is interesting in that it predicts that in 
environments in which only massive fragments are born in 
the disc, there is very little chance of finding any planet 
with mass intermediate between a few Earth masses and a 
few Jupiter masses. Fig. [T^] shows the mass of the atmo¬ 
sphere irrespectively of whether the fragment is disrupted, 
in which case the atmosphere remains bound to the core, or 
collapses, in which case the “atmosphere” simply becomes 
the part of the much more massive gaseous envelope of the 
planet. 


10 RADIAL DISTRIBUTION OF PLANETS 

We finally discuss the distribution of the simulated plan¬ 
ets over the planet - host separations. This is one of the 
most interesting observables. Unfortunately it is also one 
of the most uncertain results of our model at this point. 
While many of the planet properties studied in previous sec¬ 
tions, including the metallicity correlations, depend mainly 
on the planet formation processes of TD (fragment forma- 
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Figure 19. Mass of the core’s atmosphere versus the fragment’s 
initial mass. Note that all fragments, whether disrupted or not, 
are plotted. For non-disrupted fragments the atmosphere simply 
becomes the part of the much larger gaseous envelope and is not 
distinguishable from it. 


tion, contraction, core growth, tidal disruption, etc.), the 
radial distribution of planets depends at least as sensitively 
on the protoplanetary disc evolution model accepted here. 
Section FlU.41 shows why this statement is true for our simu¬ 
lated models, but it is also bound to be true for any planet 
formation theory fe.g.. [Alexander fc Armitagell200^ . These 
uncertainties go beyond planet formation theory and can¬ 
not be resolved until evolution of protoplanetary discs is 
understood in sufficient detail. We believe that real proto¬ 
planetary discs are far more complex than the simple model 
accepted here or in any other pub lis hed p o pulation synthe¬ 
sis models ( e.g., Ilda fc Lhil l2004al f9. l200^ : iMordasini et al.l 
l2009al . I 2 OI 2 I 1 . This view is motivated by the increasingly 
popular ideas of episodic accretion on young stars (e.g., 
iDunham fc Vorobvoy 2012h that may have direct connec¬ 
tions to GI model ( Bolev et al.ll201^ : iNavakshin fc Lodatol 
l2012li . and observational surprises in the population of 
’’transition discs” that were expected to be the link between 
protoplanetary and debris discs (see the end of the Discus¬ 
sion section). 

With this caveat in mind, fig. 1201 shows how the planet 
final separations are distributed for four different groups of 
planets. The top panel (a) shows “Earths”, planets with 
mass 0.3M 0 < Mp < 2M0. Panel (b) shows Super Earths 
defined to be planets with mass 2 M 0 < Mp < 15 M 0 . The 
third and the fourth panels show “giants” with mass between 
50 M 0 and 5Mj. Most massive gas giants, Mp > 5Mj, are 
shown in panel (d). The histograms are further broken into 
the metal rich (red), the metal poor (green shaded) and the 
total populations (blue colour). We do not show in fig. l20l the 
“assimilated” group of planets, that is those that migrated 
all the way to a = Rin, since we cannot constrain their 
further fate here. 







Figure 20. The histogram of planet final separations for four 
groups of planets: Earths, Super Earths, Giants and Super Giants, 
from top to bottom. The colours show all, metal rich and metal 
poor populations, as indicated in the legend on the top panel. 


10.1 Low mass cores 

Panel (a) of fig. [20] shows that Earth-like core-dominated 
planets are most abundant at a ~ a few AU separations. 
Physically, this is because fragments migrating from the 
outer disc are disrupted most frequently in this region (this 
can be d educed already fr om simple analytical estimates, 
see, e.g., iNavakshinIl2010i] l . and hence this is where these 
cores are deposited by the disruptions. Low mass cores do 
not migrate inward significantly, so they form the dominant 
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peak in the same spatial region. One also observes that the 
histogram for the low mass rocky planets is mainly com¬ 
posed of cores born at low metallicity environments, which 
is to be expected given the results of 31 

One may also note that the metal rich population of 
low mass cores forms a flatter distribution of host - planet 
separations which is shifted a little closer to the star. This 
can be attributed to two effects. Firstly, the cores born in 
metal rich environments are born in more compact gas frag¬ 
ments. These fragments are disrupted closer to the host star. 
Secondly, cores around metal rich stars are more massive on 
average that those in the low metallicity discs, and so they 
also migrate further in. 

10.2 Super Earths 

Panel (b) in fig. [20] indicates that super Earth planets are 
shifted towards smaller separations with respect to the low 
mass cores. This is due to a higher inward migration speed 
of the super Earths which brings most of them into the inner 
region. For this reason the super Earths dominate the mass 
function at small separations, as was already seen in fig. (6] 
The metal rich and the metal poor populations of the super 
Earths are quite similar, which again shows that there is no 
met allicity correlation f or these planets. 

ISilburt et al.l (I 2 OI 5 I I reconstructed the frequency of oc¬ 
currence of Kepler detected planets in the size range 1 —4i?0 
in the orbital period range from 10 to 200 days. The 
dashed purple line in panel (b) in £g. [20] shows their re¬ 
sult, which should approximately map to the mass range 
we use for this panel. We see that there is a discrepancy of 
about a factor of three between the observations and our one 
fragment simulations. Another way of putting this is to say 
that if different fragments did not interfere with each other 
(which is most likely not true), then we would need ~ 3 
fragments per star to reproduce the observed abundance of 
super Earths inside the inner ~ 1 AU. 

We note that the outer region of ~ tens of AU contains 
very few cores in our simulations (see also Eg. [2J ■ This out¬ 
come stems from tidal disruptions of fragments becoming 
possible only in the inner 10 AU region. This result does 
depend on assumptions and parameters of the model, such 
as the fragment migration speed, pebble abundance in the 
disc, the initial disc mass, etc. For example, simulations in 
paper II did produce massive cores in the outer tens of AU 
region (cf. fig. 4 in that paper). 

10.3 Gas giants 

Panels (c) and (d) of fig.[2n]show the same simulated popu¬ 
lation of planets with mass between 50 M 0 and 5Mj, ex¬ 
cept panel (d) is a zoom-in on the inner 10 AU part of 
the histogram. These panels show that very few of the 
gas giant planets actually end up in the region between 
0.1 < a ^ 5 AU. The most common outcome (cf. also fig. 
[SJj) for the initial fragments is to migrate inward rapidly and 
be destroyed or to collapse but end up on the inner edge of 
the disc at a = 0.09 AU. The next most common outcome is 
that some fragments are ’’stranded” in the outer region by 
the disc disappearing more rapidly than the fragments can 
migrate. 


The fraction of the simulated gas giants in the in- 
ner disc appear s too low compared with the observations. 
[Gumming et al.l (I 2 OO 8 II show that the occurrence rate of gas 
giant planets with period less than 5.5 years and in the mass 
rate 0.3 Mj to 10 Mj is « 10%. Th e frequency of t he ob¬ 
served giants around FGK stars from iGumming et al.l (l2008f) 
is shown in panel (c) of [20] with the purple dashed line. 

iBiller et al.l (I 2 OI 3 I I places a model-dependent upper 
limit of ~ 10% on the frequency of 1 — 20 Mj directly imaged 
mass planets at separations 10 — 150 AU. This shows that 
systems such as HR 8799 are very rare. Assum ing a flat in 
log a d istribution for such planets, we plot the iBiller et al.l 
(I 2 OI 3 II result as a black dotted line (recalling that we have 
20,000 host stars in our simulations). 

Panel (c) shows clearly that while our model is 
marginally consistent with the upper limits on the directly 
imaged planets, the number of planets in the inner ~ 10 
AU is very much smaller than observed. Physically, this is 
due to a very rapid inward migration of the giant planets 
through this region on the way to a = Rin, due to which it 
is quite unlikely that a planet would stop there at the time 
when the disc is dissipated away. 

The zoom-in on the inner part of panel (c) is shown 
in panel (d). It shows that the fraction of medium mass 
giant planets in our sample (which is about 1.0 %) is almost 
exactly 10 times lower than the observed gas giant fraction 
(the purple dashed curve in panel [d] is same as in panel [c] 
but is multiplied by the factor of 0.1). 

Evidently, the simulations could be reconciled with ob¬ 
servations by either a different disc model which would 
have the planets migrate slower through the inner region, 
or by having ~ 10 fragments per host star initially and 
removing the “excess” far away giants by some process. 
Glose planet-planet scat tering during the formation phase 
llVorobvov &L Basrj] I 2 OI 3 II and stel lar interactions ma y re¬ 
move some of the far-out giants llDavies et al ][2^, but 
it seems doubtful that these processes would be effective in 
removing as much as 90% of the original population. 

One feature in the radial distribution of gas gi¬ 
ant exopla nets not reproduced by power-law fits such 
as that by [Gumming et al.l (]2008|) is t he sh arp “pile-up” 
of planets at a ^ 1 AU IWright et al.l (l2009ll . Models by 
[Alexander fc Pascuccil (I 2 OI 2 II show that photo-evaporation 
of the disc may produce a similarly strong pile-up of plan¬ 
ets at about the right place. Physically, when the disc is 
dissipated by photo-evaporation, the gas is removed most 
rapidly from the region with R ^ a, few AU, which then 
implies that the planets migrate through this region slower 
than one would expect from a non photo-evaporating disc. 
This then forms a spike on the final semi-major distribution 
of gas giants. Note that this physics is independent of how 
the gas giant planets are formed, although there is a strong 
dependence on the timing of the planet “injection” into the 
inner few AU of the disc. 

Although our models include disc photo-evaporation, 
we do not find such a strong pile-up of gas giants in our cal¬ 
culations. However, further tests, to be presented elsewhere, 
show that pile-up formation is a strong function of the disc 
viscosity parameter ao. In this paper we kept it at a fixed 
“reasonable” value qq = 0.005 (cf. §2.2.1). Simulations with 
a broad range of ao actually do result in a cliff-like decrease 
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of gas giants inward of a ~ 1 AU. We plan to analyse and 
present this issue in future work. 

Finally, panel (e) of fig. [20]shows the radial distribution 
of planets more massive than 5 Mj. As for less massive 
planets, most of these are located at large separation. There 
is no strong metallicity correlation either in the inner or the 
outer disc for these planets. 


10.4 Disc evolution impact on planet yields 

As can be gathered from >12.2.H and H2.2.41 even with a fixed 
radial structure of our initial discs and a small - only a factor 
of two - range in the initial total disc mass, the time scale 
for the protoplanetary disc dissipation varies widely between 
the runs in line with the observed range in the disc lifetimes 
(fig. 1). It may be expected that the duration of the disc life¬ 
time has a significant impact on how the fragments evolve. 
To explore this issue, we define an ’’evaporation time”, which 
is an estimate of the time scale on which the protoplanetary 
disc would be removed by photo-evaporation. 


where Ma and Mev are the initial disc mass and the to¬ 
tal photo-evaporation rate, respectively. We then define the 
planet yield as the frequency of a given type of planet for¬ 
mation for four planet sub-samples: hot giants at the inner 
disc edge, e.g., all planets more massive than 50 M 0 with the 
final planet-host separation a — Rin', hot and cold giants - 
same mass range but 7?in < a < 5 AU and a > 10 AU, 
respectively; and hot sub-giants, Rin < a < 5 AU and 
Mp < 25 Mq. The exact dividing masses are not important 
for what follows. 

Figure [21] shows how the planet yields vary with the 
evaporation time. The meaning of the trends are physically 
clear. The faster the disc is evaporated, the smaller is the 
chance of a gas giant planet being ’’stranded” in the outer 
disc. In the limit of no photo-evaporation, tev —^ 00, our 
discs always manage to push the gas giants into the inner 
disc. Discs that are removed very quickly (short lev), on the 
other hand, disappear sooner than they are able to push the 
planet close to the star. For shortest tev bins, most of the 
initial fragments end up as cold giants. 

For the hot giant planets the trend is not monotonic. At 
short tev, the longer the disc is present, the more of the cold 
giants end up closer to the star, hence the yield of the hot 
gas giants increases as tev increases. However, this upward 
trend levels off and then turns into a decreasing one at longer 
tev since then the hot giants are more likely to be pushed 
all the way to a = Rin, which is clear from the dashed blue 
curve, or be disrupted and become a sub-giant planet (the 
green dot-dash curve). 

Figure [21] demonstrates that the outcome of planet for¬ 
mation process in the context of TD hypothesis for planet 
formation strongly depends on the disc model accepted, and 
this dependence is quite varied and may be non monotonic 
for different types of planets. While more experiments with 
different disc models are needed, it appears likely to us that 
the fraction of the cold giants formed in our models may be 
reduced by reducing the maximum photo-evaporation rate 
in our models (which is indeed very high, e.g., max[Mev] = 


Disc dissipation and planet yields 
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Figure 21. The frequency of cold and hot giant and hot sub¬ 
giant planet formation in our full simulated planet sample versus 
the estimated disc photo-evaporation time. The longer the disc 
exists, the smaller the chance of a cold giant planet formation, 
and the higher the chance that either a super-earth planet forms 
or that the gas giant is pushed all the way into the star (a = Ain). 

3 X 10”^ M 0 yr“^). We therefore do not consider the mis¬ 
match in the radial distribution of our simulated planets as 
an unsurmountable challenge to the TD model. 


11 DISCUSSION 

This paper presented a synthesis of 20,000 planet forma¬ 
tion experiments in the context of a relatively young planet 
formation theory called Tidal Downsizing. TD planet for¬ 
mation process begins as in GI model, with formation of a 
massive gas fragment in the outer self-gravitating gas disc. 
Instead of stopping there artificially, TD continues with (a) 
migration of the fragment in, as found in a dozen indepen¬ 
dent numerical simulations; (b) grain sedimentation forming 
the massive solid core in the centre of the fragment; (c) ei¬ 
ther collapse of the fragment, which results in a very young 
gas giant planet, or tidal disruption of the fragment once 
it migrated too close to the host star, which forms a core- 
dominated planet. TD model is thus GI plus modern physics, 
and is an attempt to rectify a certain injustice that GI model 
suffered while attention of the modellers in the last decades 
was focused on CA almost entirely. 

The most important result following from our calcula¬ 
tions is that many of the observed properties of the Solar 
System planets and exoplanets that were previously claimed 
to support GA theory uniquely are naturally reproduced by 
the TD theory as well. In particular, we find that 

(i) giant planets do contain massive cores (®. Their mass 
is dependent on the dust growth and sedimentation physics 
and the initial conditions for the fragments. At this point 


















we cannot rule out nearly zero core masses for some giant 
planets, especially for those more massive than a few Mj, 
since they con tract rapidly and may vap orise their grains too 
soon (cf. also iHelled fc Schubert! F 2 OO 8 I I. The most massive 
cores formed in the simulations are ~ 10 — 20 M@ (fig. [O]); 

(ii) gaseous envelopes of gas giants are strongly enriched 
in metals compared with their host stars (especially in or¬ 
ganics and water since rocks are more efficient in sediment¬ 
ing into the core), see m 

(iii) the over-metallicity of gas giants decreases as the 
planet mass increases, as observed for the Solar and some 
giant exoplanets (SHU; 

(iv) the model is able to produce copious core-dominated 
planets with mass from sub-Earth up to ~ 20 M®. These 
planets are the most frequent outcome of TD planet forma¬ 
tion hypothesis, but less so at high metallicities when giant 
planet survival becomes much more likely (cf. fig. [3)3); 

(v) a positive correlation of frequency of appearance of 
gas giant planets with metallicity of the host star for mod¬ 
erately massive giants in the inner ~ 5 AU region (cf. §fjS] 
and l5.ll) . 

(vi) Same planets do not show a metallicity preference 
at large separations (a > 5 — 10 AU) from their host stars 

(iED; 

(vii) The most massive giant planets do NOT correlate 
positively with metallicity (cf. > 15.31) . Also, extending the re¬ 
sults of f[8]into the higher mass regime, BDs and low mass 
stellar companions formed by GI of their discs (and then 
having accreted more gas from it) have no physical reason to 
be over-abundant in metals compared with their host stars; 

(viii) Core-dominated planets show a complicated pat¬ 
tern of occurrence versus the host star metallicity depen¬ 
dence. While the average core mass increases with metal¬ 
licity, the number of the cores formed in the disc does 
not (SST). The latter is directly connected with the fact 
that more initial gas fragments survive at high metallicities, 
which implies fewer tidal disruptions, and hence fewer core- 
dominated planets. Sub-Neptune mass planets thus correlate 
in mass but not in numbers with the host star metallicity. 
We emphasise that this prediction stems from the basic me¬ 
chanics of assembly of different types of planets in TD, e.g., 
more tidal disruptions equals less gas giants but more core- 
dominated planets. It is hence very robust, unlike the CA 
result that depends on deta il of type I migration prescrip¬ 
tions (iMordasini et al.]|2012l l: 

(ix) the composition of simulated cores is dominated by 
rock and Fe, not ices, especially at high core mass end (dTj). 
This prediction is significantly different from CA model. To 
remind the reader, core compositions are currently far from 
a settled issue even in the Solar System: while Uranus and 
Neptune are often referred to as ”ice giants”, there is no 
direct evidence that their bulk composition contains domi¬ 
nant amounts of ice. Modelling of Voyager and other data for 
these planets with a wide range of equations of state shows 
that compositions of these planets coul d be dominated by 
pure rock not ice fsee lHelled et al.ll201ll '): 

(x) due to a variety of fragment evolution histories, there 
is no one-to-one relation between the atmosphere mass and 
the core mass for the simulated planets. However, in general, 
the masses of the atmospheres of core-dominated planets are 
a tiny fraction of the total planet mass for cores less massive 
than Mcore < 5 M®, but become comparable to the core’s 
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mass for Mcore ~ 5 — 15 M® (>[6]). The atmospheres are also 
metal-rich ffig. 1151) : 

(xi) the planet mass function is dominated by the sub- 
Neptune mass planets, and contains a pronounced “tidal 
disruption desert” between Mp ~ 20 M® and Mp ~ 100 M® 
(cf. EXD; 

(xii) the planet mass function does not run away towards 
low mass (Mcore ~ IM®) cores, unlike CA mass function 

(cf. mn- 

(xiii) The PMF is smooth and does not divide on “rock- 
dominated” and “ice-dominated” cores 1 114.2.21) . 

(xiv) The PMF has a rollover at mass ~ 10 M® since it 
appears difficult to make more massive cores 1 114.2.31) . 

Point (vii) may have significant implications beyond 
planet formation theories. It is usually said that forma¬ 
tion scenarios of giant planets must be physically differ¬ 
ent from that of brown dwarfs and low mass stellar com¬ 
panio ns (e.g., see references in §2.1 of IWinn fc Fabrvckvl 
l2014fl . mainly because the metallicity correlations of these 
objects are different. While giant planet frequency of de¬ 
tection is positively correlated with metallicity of the host, 
for BDs such a correlation does not seem to exist conclu¬ 
sively, and for the low mass stellar companions it is the low 
metallicity hosts that a re more likely to host the companion 
llRaghavan et al.ll201Clh . We find that there may be physi¬ 
cal reasons within the TD formation scenario that explain 
the divergent metallicity correlations for these different ob¬ 
jects. In particular, if formation of planets starts with the 
birth of a M/ a few Mj gas fragment in a gravitation¬ 
ally unstable disc, and it does not grow more massive by 
gas accretion, then pebble accretion favours the survival of 
close-in exoplanets in high metallicity environments. If for¬ 
mation of BDs and low mass stellar companions starts with 
the birth of a more massive seed, e.g., M/ 10 Mj, then sur¬ 
vival of these by radiative cooling is favoured at low metallic¬ 
ities (see also [Helled fc Bodenheimei[l201lh . We found this 
low-opacity mode of collapse to become important already 
at ~ 10 Mj in >15.31 and it is clear that for higher mass 
fragments radiative cooling will dominate over the pebble 
accretion collapse even more. 

If stellar mass secondaries can indeed grow from low 
mass Mf ^ 10 Mj seeds born in the early embedded phases 
of the star evolution, then there are further interesting impli¬ 
cations stemming from TD model for the formation of plan¬ 
ets in stellar binary systems. Kepler observations show that 
planets orbiting binary stars may be as co mmon as ~ 1 — 10% 
JWelsh et aI.ll20f^ : ISchwamb et al.ll2013l l. This is surprising 
in the context of CA. Planetesimal velocity dispersions are 
whipped by the secondar y in stellar binary systems (e.g., 
[Paardekooper et al.l l2008tl . Planetesimals are expected to 
collide in such systems at velocities fro m tens of m/s to ~ 1 
km/s and fragment rat her than grow |[Paardekooper et al.l 
I 2 OI 2 I: [Lines et al.l l2014ll . making it challenging to explain 
the observed systems. In contrast, formation of these sys¬ 
tems does not appear to strain any physical limits in TD 
but it does require the following chain of events to occur. 
In particular, suppose that several gas fragments are born 
in the outer disc, and one of them is much more massive 
than the rest and grows by gas accretion into a stellar sec¬ 
ondary. It does not obviously precludes the other low mass 
gas fragments from maturing into planets although quantita- 
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tively their evolution is probably different from the one in a 
disc around a single star. The planet-mass fragments can be 
born inside or outside of the orbit of the more massive frag¬ 
ment. Formation of either circum-primary or circum-binary 
planets can then be achieved by shrinking the separations 
of the planet and/or the secondary due to the gravitational 
torques of the disc, assuming the disc is massive enough or 
is continued to be supplied with more matter from the outer 
envelope. 

One feature of the model which is at odds with the ob¬ 
servations is the radial distribution of the giant planets, with 
the ratio of cold (planet-host separation a > 10 AU) giants 
to warm giants (a < 5 AU) being larger than observational 
limits by a factor of 10 or mor e. Our simulations also did 
not reproduce the period valley ilWright et al.l[200^ in that 


The number of far-out gas giants can be reduced if some 
of them self-destruct due to an overly luminous core, as 
in the proposed formation rout e for Ne ptune and Uranus 
by Handburv fc Williams! (Il975h , see also lNavakshin fc Chal 
(l2012ll . In addition, the radial distribution of giants is also 
very sensitive to the disc evolution model, e.g., the very 
poorly known disc viscosity parameter a, the locations of 
fragments birth, and also the very end phases of the disc 
dispersal (cf. 1)10.411 . We believe that the current paradigm 
of protoplanetary disc evolution, which is the base of our 
disc model here, is unlikely to be entirely correct. This 
paradigm, developed with CA model in mind, all but dis¬ 
counts the importance of massive gas planets and brown 
dwarfs for the disc evolution, at best adding them post 
factum as an interesting but generally unimportant per- 
turbers. The real picture is likely to be much more com- 

J ilicated as hinted by observations of “tr ansition discs” 
Andrews et al ] I 2 OIII : lOwen fc Clarkd I 2 OI 2 II , with massive 
planets not only taking the mass from the disc but also giv¬ 
ing it back w hen these objects are disrupted by tides from 
the host star (iNavakshin fc Lodatol[20l3iNavakshinllioiSll . 
This cannot possibly be added as a post factum perturba¬ 
tion to the disc evolution since the number of gas fragments 
born in the outer disc per star may be significant (e.g., 
IVorobvov fc Ba^ l200fih and the total mass in these frag¬ 
ments dur ing the overall disc lifetim e be as much as 0.05 Mq 
or larger dCha fc Navaksh'inll201lll . This view joins up well 
with the developing paradigm of episodic accretion of young 
protostars in which all stars consume a number of mas¬ 
sive clumps ship ped into the star by disc migration from ~ 
hundreds of AU ([Dunham fc V orobvovll2012 : lAudard et al.l 
l20l4IVorobvov fc BasiJl201!t . 


12 CONCLUSIONS 

Here we presented the third in the series paper on the pop¬ 
ulation synthesis of planet formation experiments in the 
framework of the Tidal Downsizing (TD) model. A number 
of observed facts and correlations for exoplanet populations 
are reproduced by our model. On this basis we believe that 
TD model is a physically attractive alternative to CA in 
explaining many if not all of the observed exoplanets. We 
presented predictions for future observations that may dis¬ 
tinguish this theory from CA. We hope that this work will 


stimulate more theorists to contribute to the development 
of TD theory for planet formation. 
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